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ABSTRACT 

Surface  damage  to  optically  nonlinear  materials  has  been  studied 
both  experimentally  and  theoretically.  The  most  important  conclusion 
reached  to  date  is  that  laser-induced  damage  to  transparent  materials 
is  characterized  by  a  probability  for  damage  at  each  incident  power 
density.  This  is  in  contrast  to  the  previously  held  view  that  there  was 

a  threshold  power  density  which  divided  damaging  levels  from  those  which 
would  do  no  damage. 

A  model  for  the  damaging  process,  based  on  the  probability  measure¬ 
ments  and  electron  avalanche  breakdown,  has  been  devised.  In  its 
simplest  form  this  model  has  successfully  explained  the  most  imoortant 
properties  of  the  measured  damage  probability  in  over  10  different  materials. 
This  model  is  being  refined  to  examine  such  areas  as  the  dependence  of 

damage  probability  on  laser  frequency  and  the  temporal  distribution  of 
breakdown  starting  times. 

The  present  program  has  also  explained  the  difference  between  damage 
produced  by  multimode  and  TEM0(,  mode  laser  beams,  shown  that  the  form 
o  damage  (whether  on  the  surfaces  or  in  the  volume)  depends  on  the  laser 
earn  geometry  and  demonstrated,  using  streak  camera  photography,  the 
upstream  movement  of  the  beam  focal  point  under  self-focusing  conditions. 

In  addition  a  refined  theoretical  analysis  of  self-focusing  with  both 
instantaneous  and  finite  response-time  mechanisms  has  been  developed. 
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I. 


INTRODUCTION 


The  essential  component  In  many  modulators  and  Q  switches,  and 
in  all  frequency  doublers  and  parametric  oscillators,  is  an  optically  nonlinear 
crystal.  Unfortunately,  the  nonlinear  materials  having  the  largest  non- 
linearities  (i.e. ,  linear  electro-optic  coefficients,  second  harmonic  generation 
(SHG)  coefficients,  etc.)  are  often  those  most  easily  damaged  by  intense  laser 
beams.  Even  the  slightest  damage,  which  can  occur  on  the  surfaces  or  in  the 
volume  of  the  material,  can  render  the  crystal  useless.  Damage  to  optically 

nonlinear  materials  is  thus  a  major  limitation  in  the  design  of  high-power 
lasers  and  optical  devices. 


Most  reports  which  describe  new  optically  nonlinear  materials  mention 
the  occurrence  of  index-of-refraction  gradients  when  l.igh-average-power  visible 
light  passes  through  the  crystal.  In  LiNbOj  this  type  of  damage  is  thought  to 
be  caused  by  the  presence  of  crystalline  imperfections,  possibly  oxygen  vacancies 
which  trap  electrons.  Visible  light  frees  these  electrons,  which  then  move  away 
rom  the  imperfection.  The  resulting  internal  electric  fields  cause  index-of- 
-efraction  gradients  throughout  the  crystal,  by  means  of  the  linear  electro-optic 
effect.  Thus  a  coherent  process  such  as  phase-matched  SHG  will  be  severely 

limited.  However,  since  heating  the  crystal  permits  the  electrons  to  return  to 
their  initial  sites,  this  type  of  damage  is  reversible. 


On  the  other  hand,  when  the  light  intensity  incident  on  the  crystal  is 
high  enough,  irreversible  damage  occurs.  Because  of  residual  absorption,  very 
high  average  power  beams  can  heat  the  material  to  the  point  where  thermal-stress- 
induced  fracture  occurs.  The  present  program,  however,  is  concerned  with 
damage  caused  by  very  high-peak-power  pulses  of  light  to  crystals  which  do  not 
absorb  at  the  optical  wavelength.  When  exposed  to  such  pulses,  the  crystal's 
surfaces  may  become  pitted  and  charred,  and/or  the  interior  may  become  cracked 
or  filled  with  bubbles.  This  catastrophic  damage  can  be  corrected  only  by  re- 
polishing  the  surfaces  or  by  replacing  an  internally  damaged  crystal. 


1 


Th(  objective  of  this  program  is  to  determine  how  and  why  permanent 
damage  is  produced  and  to  find  means  to  minimize  or  avoid  its  deleterious  effects 
on  laser  systems.  The  work  to  date  can  be  grouped  into  three  general  areas: 

1.  Studies  of  the  residual  damage  after  irradiation  using  optical  and  electron 
microscopy;  2.  Studies  of  the  dynamics  of  the  damage  process  using  high-speed 
streak  photography;  and  3.  Determination  of  the  probability  that  a  particular  laser 
pulse  will  damage  a  particular  material.  The  following  list  summarizes  the 
major  results  reached  to  date: 

1.  Residual  damage  was  found  only  after  the  laser  irradiation  caused 
a  visible  spark  when  passing  through  the  sample. 

2.  The  importance  of  uniform  ("hot  spot"  free)  illumination  in  increasing 
the  ability  of  a  material  to  withstand  laser  irradiation  was  demonstrated  by  comparing 
TEMqq  and  multimode  damage. 

3.  Contrary  to  previously  held  opinions,  when  LiNbOg  and  KDP  are  used 
as  Pockel’  s  cell  or  frequency  doubler  crystals,  interne!  damage  is  produced  at 
the  same  or  lower  levels  of  irradiation  than  those  required  to  produce  surface 
damage. 

4.  The  laser  beam' s  focusing  conditions  determine  the  form  of  the 
residual  damage;  that  is,  whether  it  is  internal,  entrance-  or  exit-face  damage. 

5.  Internal  filamentary  damage  in  LiNbOg  was  found  to  be  composed  of 
a  series  of  very  fine,  nearly  planar  cracks  which  intersect  to  form  one  or  more 
long  lines  having  diameter  a  0.4um.  When  illuminated  with  visible  light,  diffraction 
effects  make  these  lines  appear  blurred  and  give  the  impression  of  a  "tube"  of 
damaged  material. 

6.  Experiments  designed  to  study  the  dynamics  of  damage  formation 

by  photographing  the  laser-induced  breakdown  (the  spark)  with  an  image  converter 
streak  camera  were  initiated.  The  results  of  these  measurements  showed  that 
when  light  is  focused  inside  the  material  the  first  damage  to  occur  is  near  the 
local  point  and  additional  damage  occurs  upstream  at  later  times.  All  laser- 
induced  internal  breakdowns  were  initiated  during  intervals  of  time  when  the 
level  of  irradiation  increased  or  was  maximum,  and  never  when  the  intensity 
decreased. 
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7.  A  model  for  the  results  in  6  was  developed,  in  which  self-focusing 
was  identified  as  likely  to  be  responsible  for  the  upstream  movement  of  the 
sparks  and  a  fast-response  self-focusing  process  was  shown  to  be  able  to  explain 
the  fact  that  internal  breakdowns  were  onTy  produced  when  the  pulse  intensity 
increased.  The  importance  of  beam  geometry  and  non-self-focusing  configura¬ 
tions  were  investigated. 

8.  There  is,  in  general,  no  sharp  threshold  for  laser-induced  damage. 
There  is,  however,  some  probability  that  a  particular  pulse  will  damage  a 
particular  material.  At  very  high  levels  of  irradiation  this  probability  is  unity; 
at  lower  levels  it  is  less  than  one,  and  at  very  low  levels  the  probability  for 
damage  may  be  vanishingly  small.  This  was  studied  with  both  Nd:YAG  and  ruby 
lasers. 

9.  Measurements  of  the  dependence  of  damage  probability  on  the  laser 
beam' s  power  density  for  ten  different  mate  ’ials  have  led  to  a  model  for  the 
damage  mechanism  based  on  avalanche  breakdown.  Good  qualitative  agreement 
has  been  obtained. 

10.  Using  the  streak  camera,  preliminary  measurements  of  the 
distribution  of  breakdown  starting  times  as  a  function  of  damage  probability 
have  been  obtained.  These  data  show  that  the  most  likely  time  for  a  breakdown 

l 

to  occur  is  before  the  peak  of  the  laser  pulse.  This  data  may  be  explained  using 
the  probabilistic  interpretation  of  laser-induced  damage. 
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II. 


experiments 


vary  widely  C  °Pera“°nal  charac'eristics  of  all  lasers  and  laser  systems 
y  dely,  ,s  necessary  to  describe  the  light  pulse  responsible  for  the 

between  l^T  “"ons  can  be  made 

P  ;  tie  :  ,  I  ra riaIS  31,11  the  reSUltS  °f  di«—  alters.  The  relevant 
properties  of  light  pulses  are  listed  in  Table  I  ThP  n„ic 
■  a  ctDie  1.  1  ne  pulse  energy  must  be 

oZ\"Z  Ter  tha"  S°me  arM‘rary  rClatiVe  Unlt  which  •»  meaningful 
only  to  the  investigator  who  define c=  it  ^  S 

•  .  ,  J  .  *  o  aeiines  it.  It  iS  adequate  to  give  the  power  densitv 

“  IT TT  T  /Cm  >  alS°  ®1VeS  thE  P“laa  duration  and  area 

oreover,  the  definition  of  pulse  duration  (FWHP)  or  the  time 

“  ST  Pr  n‘  ""  90  ■“  *  -  —  etc  )  mostbe 

and  geometry  of  thebearnTnllm  «  the“a)  ^  1,16  Pa“ern 

on  me  crystal  ,  these  too  are  essential. 

irradiated  su^“Z^;::ary  *°  S‘a‘e  °f  *ha  P^'a>  a"d  the 

aucb  as  polishing  00^  , T  ZT 

grease  can  lead  to  spurious  conrlncsirmc 

comments  about  the  cleaning  procedure  are  necessary.  addition,  such 

crystalline  properties  as  hardness  or  hydroscopic  nature  are  necessary  to  a  full 
description  of  the  sample.  necessary  to  a  full 


B*  Defining  the  Occurrence  of  Damage 

a  vi  M  In  '‘I  °Ur  6XPeriments'  rcsidual  damage  could  be  observed  only  after 

withSloOXSPar  T  Pr°dUCed  by  the  ‘"‘"action  of  the  light  and  the  sample.  Even 

be  d  t  t  I™*"  atl°n  n°  dama«e  in  the  irradiated  regions  of  the  sample  could 
be  detected  until  after  the  spark  was  seen.  Thus  we  chose  as  our  definition  of  the 

p° ::zze:i  draee  t  statement  ^  -  -•«  produced 

Zt  *  the  laser  PUlSe  thr0U*h  medium  was  accompanied  by  a  visible  spark 
ros, -irradiation  inspection  of  the  residual  damage,  however,  may  requir!  in 
spection  with  proper  illumination  and  magnification. 
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TABLE  I 


PARAMETERS  NECESSARY  TO  CHARACTERIZE 
A  PULSE  OF  LASER  LIGHT 


1.  Wavelength  and  spectral  width 

2.  Energy  in  joules 

3.  Pulse  duration 

4.  Transverse  mode  pattern 

5.  Longitudinal  mode  content 

6.  Beam' s  transverse  dimension  at  the  surface 

7.  Position  of  focal  point  if  any  or  the  beam  waist 

8„  Polarization  state 

9.  Pulse  repetition  rate 
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TABLE  II 


PARAMETERS  NECESSARY  TO  CHARACTERIZE 
A  CRYSTALLINE  SURFACE 


1.  Name  of  the  crystal 

2.  How  grown,  including  purity  considerations 

3.  Presence  of  twinning 

4.  Ferroelectric  domain  structure 

5.  Cleaved  or  polished  surfaces 

6.  Flatness  of  surfaces 

7.  Smoothness  of  surfaces 

8.  Cleanliness 

9.  Orientation  of  crystallographic  axes  and  light  beam  k  and 
E  vectors 

10.  Temperature 

11.  Atmosphere  surrounding  the  crystal 

12.  Dislocation  count 
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C .  The  Lasers  and  Beam -Handling  Optics 


Figure  1  shows  schematically  the  principal  features  of  the  laser 
damage  source.  The  experiments  in  this  program  were  performed  using  a 
pulse  pumped,  electro-optically  Q  switched,  NdrYAG,  or  ruby  laser.  The 
important  properties  of  these  devices  are  summarized  in  Table  III. 

The  combination  of  one  rotatable  and  one  fixed  polarizer  resulted  in 
a  variable  light  attenuator  which  was  highly  sensitive,  quite  reproducible, 
and  which  did  not  affect  the  laser  pulse' s  polarization,  spatial  distribution, 
or  duration.  If  the  fixed  polarizer  is  oriented  to  transmit  the  laser  polariza¬ 
tion  and  if  0  =  0°  is  the  angle  of  the  rotating  polarizer  which  gives  maximum 
transmission  through  this  attenuator  then  the  transmitted  intensity  at  any  other 
angle  of  rotation  about  the  beam  axis  is 

1(0)  =  bl  cos^  0 
o 

IQ  is  the  incident  light  intensity  and  b  is  the  fraction  transmitted  when  0=0°. 

D.  Initial  Experiments  and  Multimode  vs  TEMQ0  Mode  Studies 


The  light  from  a  multimode  Nd:YAG  laser  was  focused  through  a 
1  Ox  microscope  objective  onto  the  entrance  surface  of  several  optically 
nonlinear  materials  in  the  very  first  experiments.  During  these  experiments, 
we  cared  little  for  the  orientation  of  the  crystalline  axes  with  respect  to  the 
light  polarization,  and  used  pulses  which  varied  from  ~  10  to  30  nsec,  depending 
on  operating  conditions.  The  pulse  energy  was  varied  by  varying  the  laser  input 
energy.  We  were,  however,  careful  to  obtain  similar  surface  quality  in  all 
crystals  (at  least  \/4  flat  over  the  small  regions  irradiated),  treat  them  equally 
(except  for  those  which  were  hydroscopic),  use  a  low  pulse  repetition  frequency 
( i -  e . ,  1  pps  ),  and  to  employ  two  observers  to  obtain  two  independent  sets  of 
data.  The  damage  pattern  caused  by  a  focused  "multimode"  laser  on  a  LiNbO, 

O 

surface  is  shown  in  Fig.  2. 
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PBN-432 


100%  R 
Mirror 


KDP 

P.jck«l«  call 


Nd  YAG  or  ruby  lotar 
rad  in  aingla  tllipta 
pump  cavity 


a 


Intracavity  oppartura 
to  roatrict  oaeillotion 
to  TEMqq  mada.  A 
Mandenholl  wadga" 


Output  mirror 
dialactric  for 
Nd:  YAGjtopphira 
atalon  far  ruby 


Fig.  1  Laser  and  Variable  Attenuator  Configuration  for 
Damage  Studies. 


TABLE  III 


LASER  PARAMETERS 


Wavelength 

Nd:YAG 

1. 06  (im 

Ruby 

0.694  pm 

Energy  Multimode 

TEMQ0  Mode 

~100.  mJ 

1.8  mJ 

>  1.0  J 

1  .8  mJ 

Longitudinal  Mode 

Many 

Many 

Beam  Diameter  at  Output  Mirror 
Multimode 

TEMQ0  Mode 

5  mm 

1  mm 

6  mm 

0.  75  mm 

Polarization 

Linear 

Linear 

Pulse  Repetition  Rate 

1  pps 

1  pulse/5  sec 

Pulse  Duration  in  TEMQ()  Mode  12  nsec(FWHP)  12  nsec  (FWHP) 

(See  Fig.  30) 


Pulse  to  Pulse  Energy  Reproducibility  ±  7%  ±  10% 
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PBN-433 


O.lmm 


Fig.  2  Surface  Damage  on  LiNb03  Due  to  One  Pulse  from  a 

Q-Switched  Nd:YAG  Laser.  Average  power  density 
~250  MW/cm2;  dark  field  illumination. 


Multimode 
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The  list  of  materials  in  Table  IV,  ranked  according  to  their  ability  to 
withstand  surface  damage,  was  obtained  in  this  initial  phase  of  the  work.  Clearly 
the  table  is  incomplete,  it  represents  damage  from  a  "multimode"  beam  and  does 
not  reveal  any  of  the  functional  dependence  of  the  number  of  pulses  required  to 
damage  the  surface  on  the  power  density.  On  the  other  hand,  such  a  list  as 
this  is  indicative  of  the  differences  in  resistance  to  surface  damage  among  non¬ 
linear  materials.  In  addition,  since  the  surface  quality  of  KDP,  which  has  the 
highest  resistance  to  surface  damage,  was  not  as  good  as  that  of  LiNbOg,  this 
compilation  also  shows  that  the  differences  are  caused  by  the  material  and  not 
the  surface  finish. 

More  extensive  data  was  obtained  for  LiNbOg  and  BagNaNbgOj where 
we  measured  the  average  number  of  pulses  required  to  produce  damage  N, 
as  a  function  of  the  optical  power  density  on  the  crystal' s  surface.  The  data 
shown  in  Fig.  3  was  obtained.  In  this  figure  the  quantities  labeled  and  Pj 
are  defined  ass 

P  is  the  highest  power  density  of  a  particular  pulse  which  did  not 
<x  ^ 

damage  the  surface  of  the  crystal;  and 

Pj  is  the  lowest  power  density  which  will  always  damage  the  surface 
in  a  single  pulse. 

Note  that  by  definition 

P  <  P, 

oc  1 

and  so  for  values  of  power  density,  P  between  P^  and  Pj  the  number  of  pulses 
required  to  produce  damage  will  depend  upon  P.  As  we  shall  see  in  Secs.  II-H 
and  I,  the  most  revealing  experimentation  and  theorizing  concerns  this  intermediate 
range  of  power  densities.  At  this  point,  however,  it  suffices  to  define  these 
parameters  and  to  point  out  that  the  value  of  P^  ,  in  say  MW/  cm2,  can  be  used 
as  a  practical  measure  of  a  material's  resistance  to  surface  damage. 


In  this  stage  of  the  program  P  corresponded  to  the  power  density  which  caused 
no  damage  in  500  pulses.  This'concept  is  examined  further  inSec.  II-H  ir.  view 
of  the  probabilistic  interpretation. 
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TABLE  IV 


Several  Nonlinear  Materia  s  Listed  in  Order  of 


Decreasing  Resistance  to  Surface  Damage 
(Multimode  Illumination) 


P„(MW  , 

oo  c 
cm 


Potassium  Dihydrogen  Phosphate  (KDP) 

>  400 

Ammonium  Dihydrogen  Phosphate  (ADP) 

>  400 

7 -Diethylamino-  4-  Methylcoumarin 

400 

Coumarin 

10 

Ba2NaNb5015 

10 

LiNbOg 

6 

1000 


Fig.  3  Measured  N  =  N(P)  Using  Multimode,  Pulsed,  Q-switched 
Nd:YAG  Laser. 
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Since  the  laser  output  was  varied  by  varying  the  flashlamp  input  in 
these  early  experiments,  it  is  possible  that  the  laser  mode  varied  in  a  systematic 
manner,  resulting  in  these  particular  relationships.  However,  the  LiNbO 
damage  data  described  below,  which  was  obtained  using  the  laser -attenuator  in 
Fig.  1 ,  also  shows  a  similar  relationship  for  intermediate  values  of  P  through 
the  exact  functional  dependence  is  different. 

The  values  of  P^,  Pj  and  slope,  m  (see  Fig.  3)  depend  on  the  beam' s 
transverse  mode  or,  in  other  words,  on  the  transverse  intensity  distribution. 
Using  the  setup  sketched  in  Fig.  1  where  the  laser  was  restricted  to  the  TEM 
mode  and  the  power  density  could  be  varied  without  altering  either  the  pulse 
duration  or  the  polarization  state,  we  found  that  the  value  of  P  for  LiNbO  is 
350  MW/ cm  as  compared  to  6  MW/ cm2  for  a  "multimode"  piTlse.  *  In  this 
case,  Pj  =  1600  MW/ cm  .  The  damage  pit  on  LiNbOg  produced  by  a  single 
pulse  of  1 . 06M.m  light  in  a  TEMQ0  mode  and  having  P  =  1600  MW/ cm2  is  shown 
in  Fig.  4.  This  pit  is  ~  2.7  x  10  2  mm  in  diameter  and  ~2100  A  deep. 

Figure  5  shows  the  measured  relationship  between  N  and  P  for  LiNbO 
when  irradiated  by  a  pulsed,  Q -switched  TEMQQ  mode,  Nd:YAG  laser.  The  3 
pulse  duration  was  15  nsec  (FWHP)  and  it  was  focused  to  a  nearly  circular  spot 
of  diameter  <3x10  mm.  The  crystal  surface  was  polished  flat  to  \ /  10  and 
its  smoothness  figure  was  (0,0).  It  was  uncoated  and  cleaned  with  lens  tissue 
and  methanol  after  each  damage  spot  occurred.  During  these  experiments  the 
crystal  was  oriented  so  that  the  b  axis  was  normal  to  the  surface  and  the  light 
polarization  parallel  to  the  crystal' s  axis.  Note  the  similarity  in  the  functional 
relationships  between  N  and  P  obtained  in  the  "multimode"  and  TEM  mode 
experiments.  ^ 


Determination  of  th^  area  irradiated  proved  to  be  the  major  source  of 
error  in  the  measured  power  densities.  Several  measurements  of  this  quantity 
were  obtained  by  direct  measurement  of  the  area  burned  by  the  focused  laser 


Note  that  this  is  the  average  power  density  of  the  focused  TEMn,  mode  beam. 
At  the  center  the  peak  power  density  is  at  least  2  x  the  given  v§rue. 
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PBN-435 


0. 1  mm 


Fig.  4  Surface  Damage  on  LiNbC>3  Due  to  One  Pulse  from  a  TEMqo 

Mode,  Q-Switched  Nd:YAG  Laser.  Average  power  density 
~  1600  MW/cm2;  bright  field  illumination. 
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Measured  N  = 
Nd:YAG  Laser 


N(P)  Using  TEM00  Mode  Pulsed,  Q-Switched 


16 


beam  on  the  surface  of  several  different  materials  (stainless  steel,  anodized 
aluminum,  etc.)*  The  area  was  also  calculated  from  the  beam  diameter  incident 
on  the  objective  lens  and  the  effective  f  number.  These  various  determinations 
agreed  reasonably  well  (  ±  10  percent)  with  one  another,  but  in  order  to  be 
conservative  £ind  underestimate  the  power  density,  we  selected  the  largest  value 
of  the  area  obtained. 

The  difference  between  the  values  of  P  and  P  using  TEMQ()  mode  and 
"multimode"  beams  can  be  understood  by  considering  that  the  peak  power  density 
in  a  "multimode"  beam  (i.e.  ,  in  each  "hot  spot")  can  be  many  times  the  value 
estimated  from  the  pulse  energy,  its  duration  and  the  total  irradiated  area.  In 
each  hot  spot,  power  densities  as  high  as  in  a  TEMQ0  mode  are  possible  (see 
Fig.  2),  if  at  least  75  percent  of  the  "multimode"  energy  is  evenly  distributed 
among  the  several  hot  spots. 

All  the  preceding  measurements  were  performed  at  a  pulse -repetition 

frequency  of  1  pps,  at  room  temperature,  and  in  air.  No  noticeable  difference 

was  detected  in  P  for  LiNbO,  when  the  crystal  was  heated  to  200°C  or  when 
00  o 

the  ambient  atmosphere  was  N2»  Ar  or  vacuum. 

Since  the  surface-quality  question  inevitably  arises,  we  examined  the 
imperfection  density  per  unit  area  of  various  crystalline  surfaces  using  large 
angle  X-ray  scattering.  Even  the  smallest  crystalline  misorientation,  and  any 
cracks,  scratches,  or  other  imperfections,  show  up  as  areas  of  different 
scattering  intensity  in  this  type  of  measurement.  Figure  6  shows  these  Berg- 
Barrett  X-ra>  topographs  of  the  polished  surface  of  both  a  LiNbO^  and  a  KDP 
crystal.  This  procedure  enables  one  to  locate  surface  imperfections  and  to 
try  to  correlate  damage  pits  with  these  locations.  No  significant  correlations 
could  be  found. 

These  data  show  that  the  material,  and  not  the  3urface  treatment  alone, 
accounts  for  the  difference  in  damage  resistance.  KDP  offers  outstanding 
evidence  of  this;  it  cannot  be  polished  finer  than  X./4  and  smoother  than  (10,  20) 
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Fig.  6 


2  mm 


KDP 


Berg-Barrett.  Surface  Reflection  X-Ray  Topographs  of  Polished 
Crystalline  Surfaces. 
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but  the  resistance  to  damage  is  the  highest  yet  measured.  Improvement  in 
surface  quality  alone  then  could  improve  a  given  material's  resistance  to 
surface  damage,  but  it  will  not  eliminate  the  differences  between  materials. 

E.  Internal  Damage  Using  Optical  Microscopy 

In  the  experiments  described  n  II-D  the  beam  was  sharply  focused 
on  the  crystal1  s  surface  and  only  surlace  damage  was  noted.  However,  when 
the  focusing  condition  was  relaxed  slightly  and/ or  the  beam  was  focused  inside 
the  crystal,  internal  damage  was  easily  induced.  Filamentary  damage  and 
internal  pitting  as  well  as  surface/ volume  damage  were  observed  following 
irradiation  at  levels  which  in  the  earlier  experiments  did  not  induce  surface 
damage. 


To  improve  the  photographic  documentation  of  these  types  of  damage, 
photomicrographs  were  made  using  side-light  illumination.  In  this  manner 
damage  tracks  as  shown  in  Fig.  7  were  seen,  which  could  not  otherwise  be 
detected.  The  crystal  was  a  repolished  "hot"  LiNbO^*  sample  from  Crystal 
Technology,  Inc.  which  had  previously  been  used  as  a  frequency  doubler  and 
had  suffered  surface  damage.  The  tracks  run  from  one  end  of  the  crystal  to 
the  other,  are  very  straight,  and  contain  a  series  of  nearly  regularly  spaced 
regions  of  damage  extending  ~0.  04  mm  normal  to  the  length  of  the  track. 

Since  these  tracks  were  not  observed  in  other  types  of  inspection  and  do  not 
correspond  to  damage  purposely  induced,  it  is  possible  that  they  or  other  internal 
damages  are  induced  at  low  levels  of  irradiation  and  go  unnoticed  until  surface 
damage  occurs.  It  is  even  possible  that  such  internal  damage  might  be  the 
precursor  to  surface  damage. 

As  a  consequence  of  these  observations  several  other  surface-damaged 
LiNbOg  and  KDP  samples  used  in  the  past  as  Pockel' s  cells  or  frequency  doublers 
were  examined  for  internal  damage.  These  crystals  had  been  inspected  visually 
and  thought  to  have  suffered  surface  damage  only.  In  a  few  cases,  in  fact,  we  had 
the  surfaces  repolished  and  had  reused  the  crystal.  In  all  cases,  however,  within 

the  crystal,  in  the  volume  contained  between  the  damaged  areas  of  the  surfaces 
- 

"Hot"  LiNbCL  is  a  Li  rich  form  of  this  material  which  has  a  higher  a  axis  index 
matching  temperature  for  1.06  and  0.  53u  than  stoichiometric  LiNbO^. 
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Fig.  7 


1.0  mm 


1-aint  Filaments  Inside  LiNb03  Seen  with  the  Aid  of 
Sidelighting.  Large,  out-of-focus  areas  are  regions 
where  internal  damage  was  purposely  produced. 
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we  found  many  pits  and  bubbles  and  at  least  one  or  two  of  the  latter  ‘type  of 
tracks  {see  Fig.  8).  Elsewhere  the  interiors  of  these  crystals  are  free  of 
any  sign  of  internal  damage.  It  seems,  therefore,  that  internal  damage 
which  is  not  readily  detected  by  the  usual  inspection  procedures  does  occur 
in  LiNbOg  or  KDP.  In  fact,  since  the  use  of  each  of  these  crystals  was 
halted  at  the  first  sign  of  surface  damage,  the  internal  damage  was  induced 
at  the  same  or  lower  levels  of  irradiation.  To  realize  the  generality  of  this 
observation  it  is  important  to  recall  that  the  1.06am  beams  of  light  responsible 
for  the  damage  shown  in  Figs.  7  and  8  were  multimode,  unfocused  and  had 
pulse  durations  between  10  and  40  nsec  and  the  crystals  were  oriented  for 
both  Pockels-cell  and  frequency-doubler  operation. 

A  new  LiNbOg  sample*  was  prepared  with  all  six  faces  polished  so 
that  the  interior  of  th?  material  below  the  irradiated  surface  could  be  examined 
microscopically  with  both  back  and  side  illumination.  The  damages  shown  in 
Fig.  9  were  induced  with  the  beam  incident  along  the  crystal' s  b  axis  and 
polarized  parallel  to  the  a  axis. 

When  the  TEMQOmode  1.  06pm  laser  is  focused  on  the  crystal's  surface 
at  high  power  densities,  a  small  circular  pit  is  formed  on  the  surface  and  no 
internal  damage  is  observed.  At  somewhat  lower  power  (Fig.  9b)  densities 
this  type  of  surface  damage  is  sometimes  detected.  However,  at  ostensibly 
the  same  power  density  (Fig.  9c)  internal  filamentary  damage  and  surface/ volume 
damage  where  the  filamentary  damage  rises  to  the  surface  sometimes  occur. 

At  still,  lower  power  densities,  only  internal  filamentary  damage  is  detected. 

These  results  were  obtained  using  a  7 nsec  long  pulse  (FWHP)  as  compared  to 
the  15 nsec  pulse  duration  used  in  the  work  in  Sec.  II-D. 

The  output  energy  of  the  laser  used  in  these  experiments  varied  by 
±  10  percent  from  pulse  to  pulse  and  so  the  different  types  of  damage  in  Fig.  9 
induced  at  5.  3  GW/ cm  were  probably  induced  by  slightly  different  power  densities. 
However  this  level  of  irradiation  is  of  interest  because  it  represents  the  value 
at  which  the  dominant  damage  mechanism  changes.  In  the  case  under  discussion, 
the  surface  damage  mechanism  dominates  at  high  levels  of  irradiation  while  internal 
filamentary  damage  just  below  or  reaching  the  surface  is  dominant  at  lower  levels. 
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Fig.  8 


Direction  of  light  propogotion 


Edge  of  one  foce 


1.0  mm 


Internal  Pit  and  Filamentary  Damage  in  a  LiNbO0  Q-Switch 
Crystal.  J 
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We  found  that  this  change  in  damage  process  occurs  at  0.84  GW/  cm 
when  the  irradiated  area  was  four  times  that  in  Fig.  S.  Therefore,  the  two 
experiments  yield  the  same  value  for  the  "change-over"  power  or  energy  within 
experimental  error.  That  is,  when  a  small  area  of  LiNbOg  is  irradiated  by  a 
TEMqo  mode  1. 06um  laser  pulse  of  7  nsec  duration  and  —  0.5  mJ  energy  content 
the  damage  process  which  dominates  may  be  either  that  which  results  in  surface 
damage  or  that  which  results  in  internal  filamentary  damage. 

LiNbOg  frequency  doublers  and  Pockels  cells  undergo  surface  damage 
when  irradiated  by  Q-switched  multimode  Nd:YAG  laser  beams  of  ~  10 nsec 
duration  and  20-40  mJ  total  energy.  In  general,  the  cross  sections  of  these 
beams  are  composed  of  between  10  and  40  hot  spots  and  so  the  area  of  the  crystal 
irradiated  by  each  "hot  spot"  is  subjected  to  from  0.5  to  4mJ  of  1.06pm  light. 

This  is  in  good  agreement  with  the  value  obtained  above  for  the  level  of  irradiation 
at  which  th'  dominant  damage  mechanism  becomes  that  responsible  for  surface 
damage.  With  less  energetic  but  otherwise  similar  laser  beams  internal  damage 
which  generally  goes  unnoticed  is  induced. 

Filamentary  damage  has  been  associated  with  the  phenomenon  of  self- 

focusing  where  a  nonlinear  interaction  between  the  light  and  the  medium  leads 

2  3  4 

to  a  local  increase  in  the  index  of  refraction  and  so  focuses  the  beam.  ’  ’ 

If  the  focusing  is  sufficiently  strong,  the  beam  diameter  can  become  so  small 

that  diffraction  effects  become  large  and  a  trapped  filament  may  result.  The 

extremely  high  optical  electric  fields  in  the  irradiated  region  might  cause  elec- 

2 

trical  breakdown  and  thereby  damage  the  medium. 

When  the  laser  was  focused  on  the  entrance  face  of  the  crystal  the  depth 
to  which  filamentary  damage  penetrated  was  found  to  be  independent  of  the  pulse 
energy  (see  Fig.  9)  but  was  inversely  proportional  to  the  focal  length  of  the  lens. 
The  length  of  the  filament,  however,  did  vary  with  the  incident  energy.  Since 
the  light  entering  the  crystal  under  these  conditions  is  divergent,  it  is  reasonable 
that  a  short-focal-length  lens  produces  filamentary  damage  which  penetrates 
deeper  than  that  produced  with  a  long-focal-length  lens  using  the  same  incident 
power;  the  beam  must  propagate  further  into  the  crystal  so  that  self-focusing 
can  result  in  damaging  pov/er  densities. 
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When  the  lens  is  focused  inside  the  crystal,  one  end  of  the  filamentary 
damage  is  near  the  focal  point  and  the  filament  extends  back  towards  the  lens 
(see  Fig.  10).  The  end  near  the  focal  point  is  small  and  disappears  into  the 
undamaged  material  while  the  other  end  is  characterized  by  a  large  cracked 
region  having  overall  dimensions  =  10  x  the  filament  diameter.  Zverev  et  al. 
have  observed  similar  behavior  for  filamentary  damage  in  sapphire.  The 
overall  length  of  the  damage  does  not  increase  substantially  when  several 
pulses  are  fired  into  the  same  region.  However,  the  cracked  end  becomes 
larger  and  more  cracks  appear.  The  overall  length  of  the  filamentary  damage 
varies  with  the  pulse  energy,  power,  or  power  density.  As  demonstrated  in 
Fig.  11,  this  relationship  is  found  to  be  approximated  quite  closely  by  a  power 
law  where  damage  length  is  proportional  to  energy  or  peak  power  or  power 
density.  Note  that,  due  to  self-focusing,  it  is  impossible  to  give  the  power 
density  inside  the  medium;  therefore  the  pulse  duration,  energy  and  the  focusing 
conditions  assuming  linear  optics  are  specified  separately. 

A  crystal  of  KDP  which  had  previously  been  used  to  frequency-double 

a  Nd:YAG  laser  and  which  was  thought  to  be  surface-damaged  was  repolished 

for  use  in  this  work.  It  was  found  that  when  the  TEMQq  mode  NdrYAG  laser 

was  focused  on  the  surface,  damage  was  initially  produced  inside  the  crystal 

o 

at  a  power  density  of  ~  2.  7  GW/  cm  and  no  damage  at  all  could  be  produced 
at  ~  0.24  GW/cm2.  At  intermediate  levels  of  irradiation,  internal  damage 
sometimes  could  be  produced  in  one  pulse  and  sometimes  several  pulses  were 
required.  Upon  examination  of  this  sample  with  the  microscope,  however,  many 
cracked  and  pitted  regions  were  found  at  or  just  below  the  surface  which  did  not 
correspond  to  those  which  were  induced  on  purpose-  These  are  the  remains  of 
the  so-called  surface  damage  which  was  supposedly  removed  upon  repolishing. 
These  residual  damages  are  completely  invisible  to  the  naked  eye  and  so  generally 
pass  unnoticed.  Their  presence,  however,  is  responsible  for  the  erratic  results 
obtained  and  demonstrates  that  surface  and  volume  damage  occur  at  almost  the 
same  levels  of  irradiation  in  KDP  as  well  as  in  LiNbO^. 
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Direction  of  tight  propagation 

Energy  in  the 
7  nsec,  1.06/x  pulse 

0.55mJ  (5  shots) 

0.55mJ  (I  shot) 

0.40 mJ  (I  shot) 

0.27mJ  (I  shot) 

O.I7mJ  (I  shot) 

O.I7mJ  (I  shot) 

I  | 

Position  of  lens  focai  point 

Exit  surface  of  crystal -undamaged 


1.0  mm 


Fig.  10  Filamentary  Damage  in  LiNbC>3  When  Light  is  Focused  Inside 
the  Crystal  but  Near  the  Exit  Surface. 
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Fig.  11  Filamentary  Damage  Length  Versus  Pulse  Energy  Where  the 
1.06f*m  TEMoo  Mode  Laser  is  Focused  Inside  the  Crystal. 
Pulse  duration  =  7  nsec.  If  there  were  no  self-focusing,  the 
local  diameter  would  have  been  0.006  cm. 
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One  very  practical  outcome  of  this  work  stems  from  the  fact  that  any 
sample  which  shows  surface  damage  is  very  likely  also  damaged  internally; 
thus,  a  repolished  sample  is  not  as  good  as  one  which  has  never  been  damaged. 
This  holds  both  for  samples  used  in  damage  studies  and  as  frequency -doubler 
or  Pockels-cell  crystals. 

F*  Scanning  Electron  Microscope  Studies  of  Residual  Damage 
1.  Surface  damage 

Guiliano  et  al.  6  recommended  the  use  of  the  scanning  electron  micro¬ 
scope  in  the  study  of  laser-induced  surface  damage  to  ruby  crystals.  Following 
their  suggestion,  a  scanning  electron  microscope  (SEM)  was  used  to  examine 
laser-induced  damage  to  LiNbOg  with  more  resolution  and  magnification  than 
could  be  obtained  from  the  best  optical  microscopes. 

Two  SEM  micrographs  of  the  surface  damage  shown  in  Fig.  12  were 
taken  at  slightly  different  electron-beam  incidence  angles  and  a  stereo  view 
assembled.  In  this  manner  it  was  discovered  that  in  profile  the  damage  initially 
rises  sharply  above  the  undamaged  surrounding  area,  then  slowly  falls  down 
towards  the  center  where  a  small,  deeper  depression  is  found.  With  such  a 
profile,  a  surface  pit  should  more  properly  be  called  a  crater! 

The  outer  diameter  of  the  damage  in  Fig.  12  is  approximately  equal 
to  the  focal  diameter  used  to  irradiate  the  crystal.  In  the  surrounding  region 
many  microscopic  particles  are  found  which  seem  to  have  been  ejected  from 
the  crater.  Though  the  laser  was  operated  in  the  TEMQ0  mode  and  would  thus 
have  irradiated  a  circular  area,  the  crater  is  slightly  irregular  and  smooth. 

This  suggests  the  possibility  that  the  material  composing  the  crater  wall  was 
at  one  time  molten  and  flowed  into  this  shape.  In  Appendix  A  a  measurement  of 
the  absorption  in  LiNbOg  at  1.  06|j.m  is  described,  and  the  value  of  absorptivity  is 
found  to  be  too  small  to  account  for  melting.  An  additional  process  such  as  a 
rapid  increase  in  absorptivity  with  temperature  or  the  formation  of  a  hot  plasma 
which  in  turn  causes  the  damage  must  be  involved. 
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0.001  cm 

Electron  beam 
incidence  angle 
surface  normal 

0° 


0.0002  cm 

45°  from  left 
to  right 


ig.  12  Entrance  Surface  Damage  to  LiNbC>3  when  TEMqo  Mode  1.06pm 

Laser  is  Focused  on  the  Entrance  Surface.  Focal  dia.  =  0.006  cm; 
pulse  duration  =  7  nsec;  pulse  energy  =  0.26  mJ.  S.E.M. 
micrographs. 
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Figure  13  shows  the  surface  damage  caused  by  the  same  pulse  used 
to  achieve  the  damage  in  Fig.  12,  but  focused  to  a  d'ameter  one-half  as  large. 
Note  that  the  damage  crater  has  a  similar  profile  (rising  first  and  then  falling), 
the  surrounding  area  is  littered  with  ejected  matter  and  the  crater  diameter 
is  again  approximately  that  of  the  irradiation. 

The  damage  in  Fig.  14  was  caused  by  three  1. 5mJ  pulses  applied  at 
one-second  intervals.  Otherwise  the  irradiation  was  identical  to  that  used  in 
Fig.  13.  The  many  small  (more  microscopic)  pits  surrounding  the  crater  are 
probably  caused  by  molten  particles  ejected  from  the  interior  of  the  crater 
falling  back  onto  the  surface.  In  the  case  of  Fig.  13,  most  of  these  particles 
were  either  too  small  or  too  cool  to  melt  into  the  crystal.  Instead,  as  seen  in 
that  figure,  they  lie  on  the  surface.  In  this  case,  as  in  Figs.  12  and  13,  no 
internal  damage  can  be  detected  by  optical  means  in  the  region  under  the 
crater. 


Figure  15  shows  a  case  of  surface/ volume  damage  in  LiNbOg  obtained 
with  a  0.  I7mj  pulse  of  7  nsec  duration.  Note  the  very  different  character  of 
this  damage  when  compared  to  that  in  Figs.  12  -  14.  The  surface  layers  have 
been  physically  lifted  off  and  large  pieces  have  been  ejected.  A  filamentary 
damage  region  is  seen  at  the  center  of  the  outwardly  radiating  cracks.  When 
very  highly  magnified  and  overexposed,  the  SEM  micrograph  at  the  bottom  of 
Fig.  15  reveals  that  this  filamentary  damage  is  not  hollow. 

In  the  regions  nearest  the  filament  there  is  material  that  seems  to  have 
been  molten  and  to  have  re-solidified.  Apparently  a  melting  process  occurs  in 
this  case  as  well  as  in  Figs.  12-14,  but  is  obscured  by  the  much  larger  destruction 
hich  results  when  the  filamentary  damage  reaches  the  surface.  Since  the  thermal 

7 

expansion  coefficients  of  LiNbOg  are  very  anisotropic,  it  is  possible  that  follow¬ 
ing  internal  heating,  the  surface  layers  were  highly  strained  and  so  broke  off. 

The  SEM  micrographs  in  Fig.  16  shows  that  the  type  of  damage  observed 
on  the  exit  face  when  the  lens  is  set  to  focus  through  the  material  and  on  the  exit 
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Fig.  13 


Entrance  Face  Surface  Damage  to  LiNbC>3  when  TLM00  Mode 
1.06  (am  Laser  is  Focused  on  the  Entrance  Face.  Focal  dia. 

=  0.003  cm;  pulse  duration  -  7  nsec;  pulse  energy  -  0.26  mJ . 
S.E.M.  micrograph.  Electron  beam  incident  from  left  to 
right  at  45°  to  the  surface  normal. 
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Fig.  14  Entrance  Face  Surface  Damage  to  LiNbOg  when  TEMqq  Mode 
1.06  fim  Laser  is  Focused  on  the  Entrance  Face.  Focal  dia. 
=  0.003  cm;  pulse  duration  =  7  nsec;  pulse  energy  =  1.5  mJ; 
3  shots.  S.E.M.  micrograph.  Electron  beam  incident  from 
left  to  right  at  45°  to  the  surface  normal. 
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Fig.  16a  Entrance  Fact  Damage  to  LiNbC>3  when  TEMqq  Mode  1.06  |j.m 
Laser  is  Focused  on  the  Entrance  Face. 


Fig.  16b  Exit  Face  Damage  to  LiNbC>3  when  Same  Lens  as  Used  Above  is 
Set  to  Focus  the  Same  Pulse  as  Used  Above  on  the  Exit  Face. 

Focal  Dia.  =  0.006  cm;  pulse  duration  =  7  nsec; 
pulse  energy  =  0.26  mJ.  S.E.M.  micrographs 
with  the  electron  beam  incident  normal  to  the 
surface  in  both  a  and  b. 
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face  is  very  different  from  that  induced  when  the  light  is  focused  on  the 
entrance  face.  The  character  of  exit  face  damage,  as  seen  in  Fig.  16,  is 
the  result  of  a  damage  mechanism  which  is  complicated  by  the  beam' s 
propagation  through  the  medium:  Self-focusing  effects  must  be  included  in 
the  analyses  of  laser -induced  damage  whenever  the  light  enters  the  medium. 


The  form  of  damage  induced  by  a  ruby  laser  in  materials  which  do 
not  absorb  at  6943  or  3471  A  is  very  similar  to  that  produced  by  a  Nd  laser. 

When  sharply  focused  a  ruby  laser  pulse  produces  surface  damage  only;  when 
the  focusing  condition  is  relaxed,  the  damaged  regions  extend  into  the  volume 
of  the  sample.  Results  obtained  for  ruby  laser-induced  damage  are  described 
in  Sec.  II -H  after  the  concept  of  damage  probability  is  introduced. 

2*  Internal  filamentary  damage  in  LiNbQ0 

J 

Internal  filamentary  damage  in  LiNbOg  was  studied  in  more  detail  with 
the  scanning  electron  microscope  by  polishing  the  sample  in  order  to  expose  the 
damage.  The  sample  was  oriented  so  that  the  plane  of  the  polishing  lap  was 
parallel  to  the  long  axis  of  the  damage  and  material  was  removed  in  ~  5^m  steps. 
This  process  was  continued  until  no  damage  material  remained. 

Optical  observations  before  polishing  into  the  damage  indicated  that,  in 
addition  to  a  central  tube  of  damaged  material,  there  were  several  long,  thin 
cracks  lying  parallel  to  the  filament  and  oriented  so  that  they  could  intersect. 
When  the  edges  of  these  cracks  had  been  exposed,  the  scanning  electron  micro¬ 
scope  was  used  to  measure  their  width;  they  were  between  0.3  and  0.  5pm  wide. 
The  polishing  process  did  not  reveal  any  central  tube  of  damaged  material  in 
either  optical  or  electron  micrographs.  It  is  therefore  concluded  that  the  damage 
tube  observed  optically  before  polishing  is  the  overlapping  diffraction  pattern  of 
the  lines  of  intersection  of  the  several  cracks. 

The  preceding  conclusion  is  supported  by  the  fact  that  the  diameter 
of  the  "central  damage  tube"  when  measured  with  an  optical  microscope  was 
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found  to  depend  nearly  linearly  on  the  focal-length-to-aperture  ratio  (the  f  number) 
of  the  objective  lens  used.  This  result  can  be  understood  by  recognizing  that  the 
apparent  dimensions  of  a  small  object  (i.e.  ,  one  having  dimensions  on  the  order 
of  an  optical  wavelength)  when  viewed  through  an  optical  device  are  primarily 
determined  by  diffraction  effects.  In  particular,  it  is  easy  to  show  that  the 
apparent  diameter  of  a  line  source  of  true  diameter  D  is  closed  approximated 
by 

D  =  (1+1/M)  2\F  +  D 
app  ' 

when  viewed  in  light  of  wavelength  X  through  a  lens  having  magnification  M 
and  f  number,  F.  In  the  present  experiments  M  >  10,  X  «0.5(im,^6  and 
D  =  0.  4(xm,  so  that 


is  the  expected  result. 

G .  Dynamic  Properties  of  Laser -Induced  Damage 

The  sparks  which  can  be  seen  both  inside  and  on  the  surface  of  trans¬ 
parent  media  when  damaging  pulses  of  light  pass  through,  lie  along  the  beam 
axis  and  are  bright  enough  to  be  photographed  with  an  image  converter  streak 

Q 

camera.  Fersman  and  Khazov  report  that  the  temperature  of  an  entrance 
surface  spark  pmduced  by  a  Q-switched  ruby  laser  on  K-8  glass  is  initially  in 
excess  of  8OO0°K,  certainly  bright  enough  to  be  photographed  and  hot  enough  to 
cause  damage.  Since  residual  damage  in  solids  is  found  only  where  and  after 
sparks  were  observed,  these  sparks  are  either  coincident  with  or  the  direct 
cause  of  the  damage. 

In  this  work,  a  TWR  STL  Products  Model  ID  image  converter  camera 
with  streak  plug-in  units  7B  or  5B  was  used  to  photograph  the  temporal  develop¬ 
ment  of  the  laser-induced  sparks.  A  pulsed  Q-switched  Nd:  YAG  laser  which 
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could  produce  multimode  pulses  of  energy  up  to  75mJ  and  duration  of  ~  15nsec 
(FWHP)  was  the  damage  light  source.  The  peak  of  the  laser  pulse  was  reached 
in  7  -  R  nsec.  So  that  the  mode  pattern  and  pulse  waveform  would  be  nearly 
constant,  the  laser  was  always  pulsed  at  1  pps  and  with  the  same  total  input 
energy.  The  energy  incident  on  the  lens  used  to  focus  the  light  into  the  sample 
(the  "damage"  lens)  was  controlled  with  the  polarizer  attenuation  described  in 
Sec.  II-C.  A  small  part  of  the  laser  output  was  removed  with  a  beam  splitter 
and  frequency  doubled.  Several  components,  constituting  an  optical  path  less 
than  1  nsec  longer  than  that  of  the  1.06umbeam,  aimed  the  green  light  parallel  to 
the  camera  axis  and  focused  in  the  plane  containing  the  damage.  The  image  of 
this  light  provides  a  reference  on  the  streak  photo  as  to  the  time  when  light 
entered  the  medium.  Coupled  with  the  pulse  waveform,  this  streak  permits 
one  to  determine  which  part  of  the  pulse  caused  a  particular  spark.  The 
experiment  is  diagrammed  schematically  in  Fig.  17. 

Figure  18  shows  the  formation  of  sparks  in  and  on  the  entrance  face  of 
lucite  acrylic  plastic  using  two  different  lenses.  To  avoid  confusion  of  the  green 
light  streak  with  streaks  due  to  sparks,  the  green  light  was  spatially  offset  from 
the  line  of  the  sparks  by  the  indicated  distances.  When  this  offset  is  subtracted, 
these  daca  show  that  the  first  internal  spark  to  occur  in  Fig.  18a  begins  at  the 
same  instant  that  the  leading  edge  of  the  1 . 06  nm  pulse  arrives  at  the  lens  focus 
and  that  all  the  sparks  are  initiated  within  the  7-8  nsec  required  for  the  pulse 
intensity  to  reach  its  maximum.  No  sparks  are  initiated  at  later  times.  Similar 
results  are  obtained  in  Fig.  18b  using  a  different  focusing  lens.  These  photos 
show  conclusively  that,  when  the  light  is  focused  inside  the  medium,  the  sparks 
are  formed  sequentially  beginning  with  those  nearest  the  focal  point  and  ending 
with  those  nearer  the  lens.  This  result  is  also  obtained  in  other  transparent 
solids  and  liquids. 

When  the  light  is  focused  0. 123  cm  inside  the  plastic  sample,  an  internal 
spark  at  the  focus  is  the  only  one  to  appear  at  the  lowest  input  which  always  results 
in  a  breakdown.  As  the  pulse  energy  is  increased,  more  internal  sparks  appear 
until  at  ~  6/  this  input  an  entrance  surface  spark  is  also  produced.  The  delay 
between  the  arrival  of  light  at  the  focus  and  the  beginning  of  the  first  spark 
decreases  with  increasing  pulse  energy. 
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Fig.  17  Schematic  Diagram  of  Streak  Camera  Experiment. 

a.  Optical 

b.  Electrical 
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Fig.  18  Streak  Photographs  of  Laser-Induced  Sparks  in  Lucite  Acrylic 

Plastic.  In  (a)  a  3.5  cm  focal  length  lens  was  used;  in  (b)  a 
5  cm  focal  length  lens  was  used.  The  focus  was  set  to  produce 
both  internal  and  surface  sparks. 


Exit  surface  damage  can  be  produced  by  focusing  through  the  sample 
and  on  the  exit  face.  At  the  lowest  input  which  always  produces  breakdown, 
only  exit  surface  breakdown  is  observed  while  at  higher  energies  internal 
sparks  are  also  produced.  Figure  19  shows  these  sparks  in  acrylic  plastic. 

Figure  20a  shows  breakdown  produced  when  a  20 mj  multimode  1.06pLm 
pulse  is  focused  in  air.  This  spark  begins  *  3  nsec  after  the  first  light  arrives 
at  the  focus.  The  spark  in  Fig.  20b  is  formed  when  a  plastic  sample  is  placed 
0.25mm  downstream  from  the  focus.  It  begins  nearly  coincidentally  with  the 
arrival  of  the  1. 06fim  light,  and  its  growth  to  the  right  (i.e.,  upstream)  takes 
place  over  a  time  interval  nearly  2xthat  in  Fig.  20a.  In  Fig.  20c  sparks  are 
formed  both  within  and  on  the  entrance  face  with  the  light  is  carefully  focused 
on  this  surface.  The  difference  between  the  srarks  in  these  photos  shows  that 
the  surface  breakdown  which  leads  to  damage  is  dependent  on  the  presence  of  the 
material  and  is  not  due  to  breakdown  in  the  ambient  gas. 

In  order  to  study  the  fact  that  sparks  were  produced  only  when  the  laser- 
intensity  was  increasing  in  more  detail,  the  laser  cavity  was  lengthened  from 
50  to  185  cm  to  produce  pulses  having  structure  as  shown  in  Fig.  21c  and  total 
duration  on  the  order  of  100 nsec.  The  lengthened  laser  output  was  29 mJ. 

Figure  21a  shows  the  relative  spatial  position  of  the  green  light  and  the 
line  of  sparks  when  the  3.4  cm  focal  length  lens  was  used  to  damage  plate  glass. 
The  first  spark  to  occur,  as  shown  in  Fig.  21b,  began  at  tj,  14  nsec  after  the 
first  light  arrived  at  the  focal  point.  In  agreement  with  the  results  obtained 
with  the  smooth  pulse  above,  comparison  of  the  time  coordinate  for  the  formation 
of  sparks  in  Fig.  21b  with  the  pulse  waveform  in  Fig.  21c  shows  that  sparks  are 
formed  only  during  intervals  of  time  when  the  level  of  irradiation  increases. 

On  the  other  hand,  when  the  pulse  energy  or  intensity  is  high  but  nearly  constant, 
or  increasing  slowly,  for  example  between  t ^  and  t^,  no  new  sparks  are  formed. 

A  qualitative  model  which  explains  these  results  can  be  constructed  by 
requiring  that  self-focusing  of  the  beam  precede  the  initiation  of  an  internal  break¬ 
down,  that  the  self-focusing  process  is  very  fast,  say  with  response  time 
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Fig.  19  Breakdown  on  the  Exit  Face  of  an  Acrylic  Plastic  Sample 
The  lasei  was  focused  through  the  sample  and  on  the  exit 
face. 
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light  — Internal  sparks 

Fig.  20  Breakdown  on  the  Entrance  Face  of  a  Acrylic  Plastic  Sample. 

a.  Breakdown  in  air  with  no  sample  present. 

b.  Breakdown  near  the  surface  of  the  sample. 

The  focus  was  0.25  mm  outside  the  entrance 
face. 

c.  Breakdown  when  focused  on  the  sample' s 
entrance  face.  Note  internal  sparks  are  also 
produced. 
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Fig.  21  Streak  Photographs  of  Laser-Induced  Sparks  in  Plate  Glass. 

a.  Unstreaked  sparks  and  green  light 

b.  Streaked  sparks  and  green  light 

c.  Waveform  of  the  pulse  which  produced  the 
sparks  in  (b).  Arbitrary  vertical  scale. 
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~  lo"10  sec,  and  that  the  focus  moves  upstream  with  increasing  intensity.  A 
fast  self-focusing  process  can  follow  changes  in  the  pulse  intensity  and  give 
rise  to  the  sequences  of  focusing  conditions  sketched  in  Fig*  22  for  a  smooth 
pulse  and  in  Fig.  23  for  a  "wiggly"  pulse.  The  beam  focus  can  be  in  undamaged 
material  only  when  the  pulse  intensity  increases,  implying  that  new  breakdown 
can  occur  only  during  these  intervals. 

The  requirement  that  self-focusing  in  solids  have  a  fast  response  time 

o 

is  satisfied  in  glass  as  the  measurements  of  Duguay  et  al.  have  demonstrated. 
Theoretical  calculations  of  the  movement  of  the  beam  focus  in  Sec.  IV  agree 
with  the  qualitative  picture  given  above.  These  considerations  also  show  that 
the  form  and  extent  of  the  movement  is  dependent  on  the  initial  power  density 
and  geometrical  parameters  of  the  beam.  In  fact,  no  moving  focal  point  is 
found  for  a  collimated  1. 06pm,  TEMQ0  mode,  beam  of  -0.6mm  diameter  and 
—  700  MW/ cm2  power  density.  These  conditions  were  satisfied  in  an  experiment, 
and  the  results  are  shown  in  Fig.  24.  There  is  no  sequential  ordering  to  the 
initiation  of  sparks  which  is  what  one  would  expect  if  there  were  no  moving  focal 
point.  However,  once  again  all  the  breakdowns  are  initiated  during  the  first 
8  nsecs  of  the  pulse' s  duration  suggesting  the  possibility  that  the  light  which 
arrives  later  cannot  penetrate  the  existing  sparks  to  produce  new  breakdowns. 

Note  that  in  spite  of  the  absence  of  self-focusing  internal  damage  in 
plastic  at  700  MW/ cm2  is  easily  produced,  while  at  this  level  of  irradiation  and 
using  the  sharply  focused  beam  geometry  as  in  Sec.  H-H,  it  is  unlikely 
that  one  would  ever  see  surface  damage.  This  implies  that  it  is  easier 
to  produce  a  breakdown  within  a  material  than  on  its  surface.  Viewed  in 
terms  of  the  avalanche  mechanism,  discussed  in  Sec.  Ill,  there  are  more  chances 
for  an  avalanche  to  grow  to  breakdown  when  a  large  volume  of  material  is  irradiated. 
If  impurity  sites  are  important  in  providing  the  electron  which  starts  the  avalanche, 
then  this  statement  is  more  readily  understood.  The  sample  was  cut  from  a 
commercially  available  piece  of  4  ft.  y8  ft.  sheet  plastic  and  so  is  not  a  material 
from  which  in. purities  were  rigorously  excluded.  Thus,  studies  of  the  probability 
for  internal  damage  in  more  controlled  materials  are  planned  to  resolve  this 
question. 
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Sketch  of  the  Sequence  of  Breakdown  Formation  Assuming  a 
Fast  Response  Self-Focusing  Process  and  a  Smooth  Pulse 
Waveform. 
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Fig.  23  Sketch  of  the  Sequence  of  Breakdown  Formation  Assuming  a 

Fast  Response,  Self-Focusing  Process,  and  a  "Wiggly"  Pulse 
Waveform. 
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Fig.  24  Streak  Photo  of  Laser-Induced  Breakdown  in  a  Plastic  Sample 
Using  a  Collimated  TEMqq  Mode  Laser  Beam.  The  pulse 
waveform  was  not  completely  smooth.  Pulse  power  density 
700  MW/cm2.  K  J 
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Additional  information  about  laser-induced  breakdown  can  be  inferred 
from  the  length  of  time  during  which  visible  radiation  is  emitted  by  the  spark. 

The  radiant  energy  emitted  by  such  a  hot  body  in  general  rises  rapidly  as  the 
spark  is  formed  and  then  falls  over  a  much  longer  period  of  time  according  to 
the  cooling  law  for  the  particular  type  of  spark  and  medium.  Within  limits 
determined  by  its  spectral  sensitivity  and  the  level  of  initial  exposure,  the  streak 
camera  can  be  used  to  measure  the  time  interval  during  which  a  spark  radiates. 
This  provides  useful  information  about  the  mechanisms  available  to  dissipate 
the  energy  contained  in  the  various  sparks.  Since  the  sparks  in  liquids  can  cool 
by  expansion  as  well  as  by  radiation  and  conduction,  the  shortest  lived  sparks 
are  found  in  water  and  1 ,2-dichloroethane  where  they  typically  do  not  exceed 
70 nsec  in  duration.  On  the  other  hand,  the  same  laser  pulse  produced  internal 
sparks  in  plastic  and  glass  having  duration  between  500  and  2000  nsec.  Entrance 
surface  sparks  sometimes  lasted  as  long  but  generally  were  about  300  -  6000  nsec 
in  duration.  Exit  surface  sparks  did  not  exceed  300 nsec  duration. 

H.  Measurements  of  the  Probability  for  Laser-Induced  Surface  Damage 

In  Sec.  II-D,  Figs.  3  and  5  show  the  number  of  pulses  required  to  damage 
LiNbOg  and  BagNaNbgOjg  Upon  re-examination  of  these  data  it  was  realized 
that  the  measurements  could  be  more  readily  understood  if  one  accepted  the 
notion  that  the  number  of  pulses  required  to  damage  was  in  fact  a  measure  of  the 
probability  for  a  pulse  of  that  particular  power  density  to  induce  damage  in  the 
material.  Therefore,  a  new  series  of  experiments  was  performed  using  the 
techniques  described  in  Sec.  II-D  to  measure  the  probability  for  damage  as  a 
function  of  power  density  for  several  materials. 

The  experimental  setup  is  described  in  Secs.  II - C  and  D.  Throughout 
the  present  series  of  experiments  the  following  parameters  were  maintained 
constant: 
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Laser  and  laser  wavelength:  Nd:YAG  with  \  =  1.064  A  or 

Ruby  with  X.  =  0.  6943p.m 

Laser  mode:  TEMQ0,  linearly  polarized 

Pulse  waveform:  Nearly  smooth  and  symmetric  with  r  =  12  nsec  (FWHP) 

Pulse  repetition  rate:  1  pps  for  Nd:YAG  and  1  pulse/  5  sec  for  ruby 

Focusing  conditions:  A  lOxmicroscope  objective  was  used  to  focus  the 

the  beam  to  a  circular  spot  of  diameter,  d,  =  (/•003  cm 
for  Nd:YAG  and  ~  0.  0025  cm  for  ruby.  With  this 
beam  geometry  the  damage  was  confined  to  the  ma¬ 
terial  nearest  the  entrance  surface. 

Sample  temperature:  Room  temperature 

Ambient  atmosphere:  Room  air  except  for  the  hygroscopic  samples  which 

were  maintained  in  a  dry  nitrogen  atmosphere. 


All  the  surfaces  studied  were  polished  to  the  best  optical  finish  which 
could  be  obtained,  not  coated  and  carefully  cleaned  following  procedures  appro¬ 
priate  to  the  particular  material.  The  glass  and  plastic  samples,  however,  had 
lower-quality  inspection  finishes  but  were  kept  clean  and  free  of  dust  particles 
during  these  experiments. 

The  samples  were  irradiated  at  the  indicated  rate  until  the  spark  which 
coincides  with  the  occurrence  of  damage  was  observed.  The  observer  then  checked 
for  the  presence  of  residual  damage  with  a  lOOy microscope.  The  number  of  pulses, 
N,  required  to  damage  was  noted  and  then  the  sample  was  moved  so  that  undamaged 
material  could  be  studied.  At  each  power  density  about  25  different  measurements 
of  N  were  made  and  the  probability  for  damage  by  a  single  pulse  taken  to  be 

number  of  damages  ,  (1) 

P1  =  2N 

When  no  damage  could  be  found  after  500  pulses,  the  sample  was  moved  and 
irradiated  again.  If  this  occurred  five  times  in  succession,  the  sequence  of 
events  was  taken  to  indicate  that  <  0.0004.  As  discussed  in  more  detail  below, 
the  sequence  of  N  values  obtained  was  used  as  a  check  on  the  constanev  of  the 
experimental  parameters. 
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The  main  experimental  error  in  measuring  the  power  density  was 
incurred  in  determination  of  the  focal  spot  diameter.  The  true  value  of  the 
focal  spot  diameter  is  felt  to  be  within  ±  5  percent  of  that  used.  Other  errors 
enter  in  measuring  the  pulse  energy  and  duration,  so  that  the  power  density 
measurements  herein  are  estimated  to  be  accurate  to  ±  20  percent.  Note  that 
the  average,  not  peak,  power  density  is  computed  when  one  divides  pulse 
energy  by  r  and  (  */  4)d  .  Measurements  of  damage  probability  are  limited  at 
low  probabilities  by  one's  endurance  in  counting  pulses  which  do  no  damage. 

At  high  probabilities,  when  damage  occurs  within  1  or  2  pulses,  several  errors 
of  one  pulse  in  counting  the  pulses  required  to  cauae  damage  can  result  in  a 
substantial  change  in  the  value  of  Pj.  To  obtain  a  measure  of  the  error  in  p 
several  measurements  were  made  of  a  ?1  in  the  vicinity  of  0.01.  The  standard 
deviation  of  these  values  from  their  mean  was  ~25  percent. 

The  data  plotted  in  Fig.  25  shows  the  measured  probability  that  one 
of  our  laser  pulses  damages  the  surface  of  several  materials.  It  is  clear  from 
these  data  that  at  any  power  density  there  is  always  some  probability  that  a  single 
pulse  will  induce  damage.  Now  consider  the  probability,  p^,  that  damage  be 

produced  by  the  N  pulse.  If  Pj  is  the  constant  probability  that  a  single  pulse 
produces  damage,  then  p^  is  given  as 

pN  ■  <i-p,>n'1(p1)  ■  (2) 

This  is  simply  the  compound  probability  that  there  be  exactly  N-l 
nondamaging  pulses  followed  by  one  which  causes  damage.  If  a  large  number  of 
measurements  of  the  number  of  pulses  required  to  cause  damage  are  made, 
the  fraction  of  the  total  number  of  measurements  in  which  N  pulses  were  observed 
is  a  measure  of  pN-  If  the  measured  probability  distribution  and  that  predicted  by 
Eq.  (2)  using  the  measured  value  of  Pj  agree  then  one  can  conclude  that  the  prob¬ 
ability  for  damage  was  Pj  for  each  pulse.  This  also  means  pulses  of  light  were 
constant  throughout  the  experiment  and  that  the  irradiated  areas  of  the  material 
were  all  equivalent.  Figure  26  shows  the  results  obtained  for  fused  quartz 
irradiated  by  17.  9  GW/ cm'2  pulses.  Under  these  conditions  p  =  0.  16  and  99 
measurements  of  N  were  made  in  order  to  obtain  the  experimental  distribution. 
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Nd:YAG  laser  was  the  irradiation  source.  Since  the  lower 
abcissa  is  proportional  to  the  reciprocal  of  the  electric 
field,  this  plot  serves  to  test  for  a  dependence  of  the  form 
p  ~  exp  (-  const/E). 
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Bars  s  Measured  distribution 
— •—  -  Calculated  using  measured 
P  in  P  =P.  (I-P,)N_I  for  fused 
quartz  at  17.9  GW/cm2 
at  1.06/jl 


Probability  that  Damage  Occurs  on  the  Nth  pulse  Versus  N 
for  the  Fused  Quartz.  A  17.9  GW/cm2  TEMqo  mode  1.064pm 
beam  was  used. 


It  IS  evident  from  the  analysis  and  data  above  that  it  is  quite  possible 
to  have  one  sample  of  one  material  be  damaged  by  one  number  of  pulses  and 
another,  identical  sample  be  damaged  by  another.  In  fact,  for  power  densities 
where  Pj  is  small  (p^  <  0.01)  the  probability  for  not  causing  damage  in  N-l 
pulses  is  large,  10  therefore  pN  =  Pj  for  many  values  of  N,  and  it  is  highly 
likely  that  measurements  of  N  will  yield  many  widely  different  results.  Having 
made  measurements  of  N  at  this  level  of  irradiation  one  must  then  examine 
the  distribution  of  N  values  in  order  to  be  certain  that  Pj  remained  constant 
throughout  the  experiments. 

Consider  the  data  shown  in  Table  V.  The  LiNbOg  crystal  studied  in 
this  experiment  was  obtained  from  the  Union  Carbide  Company, 11  and  was 
the  only  sample  to  show  any  extreme  difference  between  the  measured  distribu¬ 
tion  and  that  calculated  from  p  .  Light  was  incident  along  the  c  axis  of  this 
sample.  At  3.15  GW/ cm2,  the  material  withstood  500  or  more  pulses  four 
times  in  a  set  of  25  measurements.  The  other  21  times,  damage  was  achieved 
in  the  first  few  pulses.  Since  the  data  in  Fig.  26  lends  strong  evidence  to  the 
notion  that  any  laser  pulse  was  almost  identical  with  any  other,  these  data 
suggest  that  this  sample  was  not  the  same  everywhere.  If  the  four  sets  of 
500  non-damaging  pulses  are  included,  then  one  finds  Pj  =  0.01.  Among  the 
25  measurements,  however,  there  are  too  many  low  values  of  N  and  too  many 
at  values  of  N  >  500  for  the  set  of  25  N  values  to  be  considered  a  valid  sample 
of  the  distribution  pN  =  (0.  99)N_1(0. 01).  If  the  four  measurements  of  no 
damage  after  500  or  more  pulses  are  considered  to  represent  an  unusually 
damage-resistant  material  and  the  other  21  measurements  are  considered  as 
representative  of  a  more  easily  damaged  main  component  of  the  sample  then, 
for  this  component,  Pj  =  0.35.  The  distribution  of  these  21  values  of  Nj  though 
obtained  from  a  small  number  of  samples,  is  close  to  that  predicted  by 
PN  =  (0.  65)N  ^(O.  35).  Therefore  we  conclude  that  this  sample  of  LiNbCL 
contams  regions  which  are  more  resistant  to  laser-induced  damage  than  the  rest 
of  the  material.  Initial  measurements  show  that  the  dimensions  of  these  regions 
can  be  ~  1  mm.  The  bulk  of  this  sample,  however,  has  slightly  higher  damage 
probability  at  a  particular  power  density  than  that  shown  in  Fig.  25  for  light 
incident  along  the  a  axis  of  a  LiNbOg  sample  obtained  from  Crystal  Technology,  Inc. 
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RawData  for  Surface  Damage  to  LjNbO^  as  a  Function  of  Power  Densit 
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Light  was  incident  along  the  c  axis.  The  crystal  was  prepared  by  Union  Carbide  Co. 


Table  VI  lists  two  operationally  interesting  levels  of  1.  06pm  irradiation 
for  several  materials.  One  is  P^,  the  lowest  power  density  at  whioh  the 
probability  for  damage  in  a  single  pulse  is  unity,  as  defined  in  Sec.  II-D.  The 
other  is  called  P2500  anC*  *s  a  Power  density  a*  which  500  pulses  induced  no 
damage  in  five  successive  trials.  This  sequence  of  events  strongly  suggests 
a  very  low  probability  for  damage. 


An  improved  experimental  arrangement  was  introduced  to  permit 
graphical  recording  of  both  the  laser  pulses  and  the  occurrence  of  damage. 

One  pen  of  a  dual-pen  chart  recorder  was  driven  with  a  signal  proportional 
to  the  laser  pulse  energy  and  the  other  pen  monitored  the  intensity  of  a  beam 
of  light  transmitted  through  the  region  which  was  exposed  to  the  laser  pulses. 
This  arrangement,  shown  in  Fig.  27,  was  sensitive  to  the  appearance  of  surface 
damage  craters  «  0.003  cm  in  diameter,  the  smallest  damage  produced.  A 
typical  recorTng,  derived  from  this  arrangement,  is  shown  in  Fig.  28  where 
a  ruby  laser  at  a  9 GW/ cm  was  used  to  irradiate  an  x-cut  crystalline  quartz 
sample.  Note  the  ±  10  percent  laser  pulse  energy  variation.  The  first  pulse 
to  strike  the  sample  is  marked  with  an  arrow  and  the  99th  pulse,  marker  with 
a  star,  is  the  one  which  produced  damage.  This  is  indicated  by  the  sharp  re¬ 
duction  in  the  transmitted  monitor  intensity.  The  laser  trigger  signal  was 
picked  up  by  the  photomultiplier  and  so  results  in  a  fiducial  marker  on  the 
record  of  the  transmitted  light  intensity. 


The  data  in  Fig.  28  shows  that  at  this  level  of  irradiation  the  sample 
does  not  have  to  damage  when  exposed  to  any  particular  prise.  In  fact,  preceding 
the  damaging  pulse  (No.  99),  there  were  several  more  energetic  pulses  which, 
according  to  the  "threshold  point  of  view,  "  should  have  been  the  ones  to  cause 
damage. 


In  Fig.  29  the  results  of  damage  probability  measurements  for  LiNLO,,, 
crystalline  quartz  and  KDP  using  ruby  laser  irradiation  are  shown  in  comparison 
with  data  obtained  with  the  Nd:YAG  laser.  This  comparison  reveals  that  the 
dominant  damage  mechanism  appropriate  to  LiNbO^  using  ruby  irradiation  is 


55 


TABLE  VI 


Two  Power  Densities  of  Interest:  =  Lowest  Power  Density 

at  Which  Damage  Always  Occurs  in  a  Single  Pulse 

Pteqq  =  A  Power  Density  Where  the  Probability  for 

% 

Damage  In  a  Single  Pulse  is  Less  than  0.  0004 


Material 

P1 

P2500 

Fused  Quartz 

24.  0  GW/cm2 

8.  3  GW/i 

Plastic 

16.  1 

11.0 

Glass 

14.4 

6.7 

KDP 

14.4 

2.  1 

LiNb03 

11.  1 

2.  0 

ADP 

6. 4 

2.0 

Crystal  Quartz 

6.  4 

1.  2 

Ba2NaNb5Oj5 

6.  4 

1.2 

SiTi03 

6.  4 

0.08 

LiIC>3 

3.  2 

0.  37 

Pj  for  fused  quartz  is  obtained  by  extrapolating  the  data  plotted  in  Fig.  25. 
All  other  values  are  measured.  Note  that  P.,50Q  is  not  necessarily  the 
power  density  at  which  pj=  0.  0004. 
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Fig'  27  aS  Diagram  °f  A  ,,0ma,ic  Pul-  »nd  Damage  Monitoring 
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Fig.  28  Record  of  Exposures  of  X-Cut  Crystalline  Quartz  to  Pulses 
of  Ruby  Laser  Light  at  ~  9  GW/cm2. 


(Power  density) 


Fig.  29  Comparison 

Obtained  Usi 
(See  also  Fig 


quite  different  from  that  using  a  Nd:YAG  laser.  When  using  a  ruby  laser  to 
irradiate  LiNbOg  the  material  is  less  likely  to  damage  than  when  using  a 
NdsYAG  laser  of  the  same  power  density.  However,  the  "slope"  of  the 
relationship  between  damage  probability  and  power  density  plotted  for  LiNbOg 
is  much  steeper  in  the  former  case.  In  fact,  with  so  steep  a  slope  this  material 
with  ruby  laser  irradiation  may  be  said  to  exhibit  a  "damage  threshold.  "  This 
case  requires  further  study  and  analysis  because,  first,  ruby  laser  light  can 
produce  annealable  optical  inhomogeneities  in  LiNbOg  in  the  manner  described 
in  Sec.  I  and,  second,  LiNbOg  absorbs  strongly  at  0.  347|i.m, the  ruby  laser 
second  harmonic  and  so  two-photon  absorption  may  be  important.  In  Secs.  III-B 
and  C  the  dependence  of  damage  probability  on  the  laser  frequency  is  considered. 
Both  the  rate  of  change  of  probability  with  optical  field  strength  and  the 
magnitude  of  the  probability  for  damage  at  any  field  are  dependent  on  the  laser 
frequency.  Cases  such  as  those  of  crystalline  quartz  and  KDP  will  \>e  studied 
in  more  detail  in  the  future  to  determine  how  to  account  for  the  data  in  Fig.  29. 

Since  the  ruby  laser  damage  data  was  obtained  during  the  last  month 
of  this  program,  there  has  not  been  enough  time  for  more  complete  experimental 
and  theoretical  studies  of  the  frequency  dependence  of  damage  probability.  This 
will  be  a  major  area  for  future  work. 

I.  Measurements  of  the  Distribution  of  Breakdown  Starting  Times 

The  experimental  apparatus  described  in  Fig.  17  with  the  Nd:YAG 
laser  in  F'g.  1  was  used  to  study  the  distribution  of  surface  breakdown  initiation 
times  as  a  function  of  damage  probability.  Since  these  experiments  were 
initiated  during  the  summer  of  1971  only  two  materials,  plate  glass  and  SrTiOg, 
had  been  studied  at  the  time  this  report  was  written.  So  far,  however,  we  have 
found  that  the  most  likely  time  for  a  breakdown  to  begin  is  before  the  peak  of  the 
laser  pulse  and  that  the  width  of  the  distribution  increases  with  decreasing  optical 
electric  field  strength.  These  results  are  discussed  analytically  in  terms  of  the 
probabilistic  model  of  damage  in  Sec.  III-B. 
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Figure  30  shows  the  distributions  of  breakdown  starting  times  obtained 

for  plate  glass  at  three  different  levels  of  irradiation.  Twenty  to  twenty-five 

measurements  were  made  to  obtain  each  distribution.  At  the  highest  level, 

where  the  peak  optical  field  strength  was  41.8  x  107  V/m  the  power  density 

was  Pj,  the  lowest  value  at  which  damage  always  occurred  in  one  pulse.  Notice 

the  spread  in  the  distribution  as  the  field  strength  is  lowered  and  the  fact  that 

the  most  probable  starting  time  is  before  the  field  strength  is  maximum  in  all 

three  cases.  Figure  31  shows  similar  results  for  SrTiO„. 

o 

The  spread  in  the  starting  times  as  the  field  strength  is  lowered  is 
further  evidence  of  the  probabilistic  nature  of  laser  induced  damage.  A 
phenomenon  governed  by  a  threshold  field  strength  would  always  break  down 
when  that  field  strength  was  reached  and  the  distribution  of  starting  times 
would  always  be  very  sharp. 
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Typicol  pulse  waveform,  Nd:YAG  laser 
(smoothed  slightly) 
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Breakdown  starting  tirro  ,  t  irsec) 

Fig.  30  Distribution  of  Surface  Breakdown  Starting  Times  for  Plate 
Glass  for  Three  Different  Damage  Probabilities. 


Typical  pulse  waveform,  Nd.-YAG  laser 
(smoothed  slightly) 
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HI.  THEORETICAL  EXAMINATION  OF  DAMAGE  PROBABILITY 


A.  Avalanche  Breakdown 


Avalanche  breakdown  has  been  considered  previously  by  several 

other  workers,  including  Wasserman, 12  Sviridenkov,  13  Zverev  et  al. 14 
1 5 

and  Hellwarth.  Howe^  er,  since  the  results  of  their  work  showed  that 
the  calculated  breakdown  field  was  an  order  of  magnitude  larger  than 
that  which  was  measured,  this  mechanism  is  not  generally  accepted  as 
responsible  for  laser  induced  damage.  The  data  presented  in  Figs.  25 
and  29  and  the  discussion  of  avalanche  breakdown  presented  in  the  rest 
of  this  section  show  that  avalanche  breakdown  should  be  considered  as 
a  possible  d 'image  mechanism. 


Our  review  of  previous  work  on  the  possibility  of  avalanche  break¬ 
down  in  optical  field  begins  with  a  discussion  oi  the  role  of  electron 
collisions.  A  perfectly  free  electron  in  an  optical  frequency  electric 
field  will  simply  oscillate  back  and  forth  with  its  velocity  90°  out  of 
phase  with  the  driving  field  and  consequently  there  will  be  no  average 
energy  absorption.  The  coherently  oscillating  energy  of  such  an  electron 
in  a  field  at  the  damage  threshold  is  readily  calculated  to  be  only  10”3  to 

a  ^ 

10  eV.  Therefore  there  is  no  possibility  of  an  ionizing  collision. 


However,  electrons  in  a  solid  are  not  really  free.  They  experience 
collisions  with  phonons  which  tend  to  relax  the  distribution  towards 
equilibrium.  These  collisions  also  result  in  a  component  of  electron 
velocity  in  phase  with  the  field  so  that  there  can  be  an  average  energy 
absorption.  A  ohenomenological  equation  of  motion  for  an  electron  in 
a  solid  may  be  written  as 


m 


TF 


t>  ■ 


eE  e 


iu  t 


(3) 
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where  m  is  the  effective  mass,  v  is  the  electron  velocity,  t  is  some 
average  relaxation  time,  e  is  the  electronic  charge  and  E  is  the  amplitude 
of  the  electric  field  at  radian  frequency  u.  It  is  easy  to  see  that  the 

energy  £,  of  an  electron  raised  to  the  conduction  band  increases  initially 
as 

i  *  2_2,„  *  2^,2 

a  jr/at  =  \  Re  <eE  •  V %  =  J-e-.E  /2m  *  — e_L _  .  {4) 

1+  u>  t  2m*u  t 

Equation  (4)  predicts  that  £  will  grow  without  limit,  which,  of 

course,  is  unrealistic.  The  energetic,  or  "hot",  electrons  will  also 

lose  energy,  primarily  by  collisions  with  optical  phonons.  An  equilibrium 

may  be  thus  established  as  discussed  by  Wasserman.  However,  at  high 

fields,  all  electrons  having  energy  greater  than  some  value  £  will  on 

c  ’ 

the  average,  gain  move  energy  than  they  lose  and  therefore  will  be 
accelerated  to  ionizing  energies. 

Detailed  treatments  based  on  this  approach  must  include  a  calculation 
of  the  electron-optical  phonon  interaction  and  a  discussion  of  the  choice  of 
£ Q.  The  problem  has  much  in  common  with  the  dc  breakdown  of  dielectrics, 
where  fairly  accurate  calculations  of  the  threshold  are  possible.  16  For  the 
optical  case,  this  procedure  was  carried  out  by  Wasserman.  12  He 
calculated  a  breakdown  field  an  order  of  magnitude  greater  than  is  observed. 


On  the  other  hand,  a  rather  different  approach  has  been  successful 
for  the  closely  related  problem  of  avalanche  multiplication  breakdown  in 
semiconductor  p-n  junctions.  There  Chynoweth17  found  empirically  that 
the  ionization  coefficient  a  (E),  analogous  to  the  first  Townsend  coefficient^ 
in  gas  breakdown,  varied  with  field  as 

a  (E)  ~  exp  (-const  /E)  (5) 

Physically,  a  (E)  is  the  average  number  of  hole-electron  pairs  produced  by 
an  electron  falling  1  cm  in  the  field  direction. 
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For  high  frequency  cw  fields,  breakdown  occurs  when  the  rate  of 

increase  of  the  number  of  electrons,  simply  related  to  a  (E),  exceeds 

the  rate  of  loss  of  electrons  due  to  recombination,  trapping  or  diffusion 

1 8 

out  of  the  field  region.  For  pulsed  fields,  there  is  the  additional 
requirement  that  the  electron  concentration  buildup  to  destructuve  levels 
during  the  pulse  period.  In  either  case,  the  major  problem  is  the 
calculation  of  a  (E)  which,  to  our  knowledge,  has  not  been  carried  out 
for  optical  fields. 


For  dc  fields  several  calculations  of  a  (E)  art  available.  The 

19 

simplest  is  due  to  Shockley.  He  contends  that  an  equation  of  motion 
like  Eq  .  (3)  is  not  applicable  because  it  gives  the  behavior  of  an  average 
electron.  Ionization  is  produced  mainly  by  those  exceptional  electrons 
which  are  accelerated  to  the  ionizing  energy  £ j  without  undergoing  a  single 
collision,  even  though  they  must  cover  a  distance  of  many  mean  free  paths. 
The  probability  of  a  particular  electron  covering  a  distance  x  without  a 
collision  is  exp  (-x/i  ),  where  A  is  the  mean  free  path.  To  be  accelerated 
to  £ .,  the  electron  must  traverse  a  distance  x  =  £  ./eE.  Therefore 

a  (E)  ~  exp  (-jP./eEi  )  .  (5a) 

Shockley  found  that  Eq  .  (5)  gave  a  good  fit  to  Chynoweth' s  data  with 
i  ~  100A. 


20 

A  very  different  approach  was  taken  by  Wolff  who  used  a 
Boltzmann  equation  and  expanded  the  electron  distribution  function  in 
Legendre  polynomials  retaining  only  the  zeroth  and  first-order  terms. 
He  found 


o 


~  exp  (-const/E^) 


(6) 


The  basic  difference  between  the  approaches  of  Wolff  and  Shockley  is 
that  the  former  assumes  that  the  distribution  is  only  slightly  distorted, 
while  the  latter  assumes  that  it  develops  a  large  spike  in  the  field 
direction. 
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21  22 

The  discrepancy  was  resolved  by  Baraff  and  Keldysh  who 
showed  that  Eq.  (5)  is  the  low-field  limiting  form  of  the  general  solution, 
while  Eq.  (6)  applies  only  at  high  fields.  The  crossover  point  isn't  well 

g 

specified,  but  Shockley' s  result  seems  to  work  well  to  E  ~  10  V/cm. 

In  other  words,  the  electron  may  make  several  collisions  before  ionization 
without  a  (E)  showing  a  significant  departure  from  the  exp  (-const/E) 
behavior  which  was  derived  on  the  assumption  of  no  collisions.  Some 
justification  for  this  success  of  Shockley' s  theory  has  been  given  by 
Baraff.  ^  * 

Shockley' s  aigument  can  also  be  applied  to  high  frequency  fields. 

For  an  electron  to  continually  gain  energy  from  an  alternating  field,  it 
must  undergo  certain  favorable  elastic  collisions  which  reverse  its 
momentum  when  the  field  reverses.  Then  a  lucky  electron  in  an  optical 
frequency  field  is  one  which  undergoes  just  the  right  collisions  to  keep 
it  in  phase  with  the  field.  Such  an  electron  requires;  only  about  5  cycles 
of  the  field  (at  E  =  Ejj  see  Table  VIII)  to  reach  e  ..  Admittedly,  the 
sequence  of  special  collisions  has  a  low  probability  of  occurrence,  but 
so  does  the  collisionless  sequence  postulated  by  Shockley.  In  fact,  if 
the  probability  of  some  collision  during  one  half-cycle  of  the  field  is 
almost  unity,  then  no  collision  (or  one  with  only  a  small  momentum  transfer) 
is  as  specific  an  event  as  the  particular  collision  needed  to  keep  the  electron 
in  phase  with  the  field.  We  would  therefore  expect  the  dc  and  optical 
ionization  probabilities  to  show  the  same  field  dependence. 

To  make  this  argument  more  precise,  suppose  that  an  electron, 
starting  from  rest,  will  remain"  acceptably"  in  phase  with  the  field 
provided  it  travels  a  distance  Xj  -  Ax1  during  one  half-cycle,  undergoes 
a  "proper"  collision  somewhere  between  x^  -  Ax^  and  x^  +  Axj,  then 


Note  that  by  using  the  Shockley  model  for  breakdown  we  are  considering 
the  energy  distribution  of  electrons  to  be  distorted  with  a  large  spike  at  high 
energies  and  in  the  direction  of  the  field.  This  is  the  same  conclusion 
reached  by  Giuliano  et  al.  using  several  physically  reasonable  approximations 
to  analyze  the  kinetic  equation  for  electrons  interacting  with  phonons  in  an 
intense  optical  field.  23 
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The  initial  free  electron  concentration  calculated  from  Eq.  (7)  setting  p(Ej)  =  1  and  assuming  f  = 
Note  that  n  =  N  4-  focal  volume. 


travels  to  x2  -  Ax2  without  collision,  etc.  Obviously  the  turning 

points,  x  ,  must  be  given  by  the  positions  of  the  electron  at  the 

times  when  the  field  reverses.  Suppose  further  that  of  all  the 

electrons  which  undergo  some  collision  between  x^  -  Axm  an<^ 

x  +  Ax  ,  a  fraction  f  undergo  a  collision  which  reverses  their 
m  m  m 

momentum  to  a  sufficient  degree.  For  isotropic  scattering,  fm  might 

be  expected  to  be  about  0.  1  or  0.  2.  The  overall  probability  per  unit 

time  of  the  electron  being  accelerated  to  €•  is  the  trial  rate,  or 

1  1 

frequency  of  collisions,  r  jj"  ,  times  the  probability  of  this  sequence 
of  events: 


P(e.) 


-1 

Tcoll 


M 

2  [exp  (- 
m=  1 


x 


m-1 


l/Ol 


[l-exp(2Axm/i)]  fm. 


(7) 


M  is  the  approximate  number  of  half  cycles  of  the  field  required  for  the  elec¬ 
tron  energy  to  reach  e  Now  if  £  is  sufficiently  small,  exp  (-2Axm/i  )<<  1, 
sc  that  any  electron  which  makes  it  into  the  range  xm±  Axm  undergoes  some 
collision  with  high  probability,  then 


M 


p<«i>-  W1  I  [exp-(|xm-xm-ll/i)]  fnT  Tcoll'  f  exP(‘VJ) 

m- 1 

where  x .  is  the  total  path  length  needed  to  attain  energy  e  .,  i.e. , 

M 

xi  s  2  lxm  ^m-l1  =  ci/eE 


(8) 


(9) 


m=  1 


with  TT  being  the  time  average  of  the  absolute  magnitude  of  the  optical  electric 
field.  We  have  assumed  that  all  f  =  f.  This  treatment  is  clearly  over¬ 
simplified  because  other  sequences  of  collisions  which  include  some  un¬ 
favorable  events  may  also  lead  to  ionization.  In  fact,  its  inadequacies  are 
the  same  as  in  Shockley' s  dc  model,  which  nonetheless  has  been  successful 
in  describing  the  experimental  data. 
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It  still  remains  to  connect  P(e.>,  the  probability  of  a  single  electron 
re.  ching  e  with  p^E),  the  probability  of  damage  in  a  single  pulse.  In 
general,  this  is  a  difficult  problem,  requiring  the  calculation  of  the 
distribution  of  avalanche  sizes  and  integrating  to  determine  what  fraction 
of  the  avalanches  exceed  the  critical  size  for  damage.  However,  the 
experimental  data  suggests  that  Pj(E)  oc  P(e.),  which  could  happen  if  the 
avalanche  statistics  were  governed  by  the  first  stage  or  two  of  the 
avalanche.  This  assumption  is  plausible  since  the  mean  energy  of  the 
secondary  carriers  produced  by  an  ionizing  collision  is  much  greater 
than  kT  and  may  even  be  greater  than  the  bandgap.  24  The  secondary 
then  has  a  much  higher  chance  of  itself  producing  ionization  than  a  thermal 
carrier,  and  so,  following  the  first  one  or  two  steps  in  the  avalanche, 
the  process  can  continue  until  breakdown  occurs.  25  Thus,  assuming 
further  that  the  mean  number,  N,  of  initial  free  electrons  in  the  focal 
volume  is  independent  of  E  prior  to  onset  of  the  avalanche,  we  obtain 


Pi  (E)  = 


..M  . 

t  N  (t 


/t  ) 
laser'  coll 


.-e./eEi 


(10) 


Now,  M  is  given  by 

n/  2me  . 

M  = - i. 

ell 


(11) 


where  vQpt  is  the  optical  frequency.  Then  the  fM  factor  may  be  written 


(12) 


so  that  this  factor  simply  changes  the  coefficient  of  T?  ~ 1  in  the  exponent. 
Therefore,  finally, 


Pj(E)  =  A  e'K^ 


(13) 


where 


A 


N  T 


laser 
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(14) 


and 


K  = 


e  . 
1 

e  i 


V2me 


2  v  ,  A  n  T 
opt  f 


The  second  term  in  K  was  omitted  in  previous  reports,  but  can  be 
significant  in  materials,  such  as  SrTiO^,  which  have  small  measured 
values  of  K. 


(15) 


The  data  in  Figs.  25  and  29  and  Table  VTI  agree  qualitatively  with 
Eq.  (13)  in  several  respects.  Obviously  the  functional  dependence  on 
field  is  correct.  Also,  the  amorphous  materials  are  all  quite  damage- 
resistant  and  show  large  slopes.  In  terms  of  Eq.  (10),  this  is  probably 
due  to  the  small  value  of  A  ,  which  should  be  comparable  to  the  mean 
interatomic  distance,  as  observed  (see  Table  I).  Among  the  crystalline 
materials,  those  with  a  large  gap  (KDP,  ADP  and  quartz)  and  therefore 
large  e  ^  have  a  large  slope.  SrTiO^  is  known  to  have  a  relatively  large 
mobility  which  implies  a  large  A  (to  the  extent  that  "hot"  electrons 
are  scattered  by  the  same  mechanism  as  thermal  ones)  and  therefore  a 
small  slope.  The  values  of  A  in  Table  VII  derived  neglecting  the  second 
term  in  Eq.  (15),  for  all  of  the  crystalline  materials  are  in  the  range 
20  -  220A.  For  comparison  A  ~  50  -  100A  is  found  in  Si  and  Ge  from 
avalanche  multiplication  data.  16,17 

A  further  check  on  Eq.  (10)  is  obtained  by  computing  N  using 
1  ancl  assuming  f  =  0.1.  The  resulting  carrier  concentrations, 
n,  shown  in  the  last  column  of  Table  VII  are  all  reasonable.  Only  the 
values  for  LiNbOg  and  plastic  look  at  all  unusual  for  dielectrics,  and 
in  the  latter  case  especially,  n  is  very  sensitive  to  the  precise  value  of 
a  10  percent  change  in  E^,  comparable  to  experimental  error,  leads 
to  a  factor  of  10  change  in  n. 

In  spite  of  this  good  qualitative  agreement,  the  model  which  led  to 
Fq.  (10)  is  highly  simplified  and  further  theoretical  work  is  needed.  Tiie 
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factors  f  should  be  calculated  taking  account  of  band  structure  and 
phonon  spectra.  The  variation  of  electron  mean  free  path  and  effective 
mass  with  energy  should  be  considered,.  Finally,  a  proper  calculation 
of  P(e  j),  possibly  by  a  Monte  Carlo  technique,  should  be  performed. 
This  would  obviate  the  need  for  the  assumption  of  small  i  which  was 
made  in  our  derivation,  and  which  in  fact  is  not  satisfied  for  some  of 
the  materials  studied  here. 

B.  Temporal  Effects 


The  discussion  in  the  last  section  considered  the  optical  field  to 
be  constant  over  the  pulse  width  Tjaser*  This  assumption  must  be 
removed,  and  a  more  rigorous  derivation  of  the  damage  probability 
given,  in  order  to  explain  the  streak  camera  data  in  Sec.  II-I. 


From  the  discussion  in  the  last  section,  the  probability  per  unit 
time  of  any  one  electron  being  accelerated  to  e  .  is  given  by  Bexp(-K/E) 
where 


B  = 


(16) 


Therefore  the  probability  per  unit  time  of  some  electron  reaching  e .  is 

h(t)  =  n^  B  dv  (17) 

where  n  is  the  density  of  free  carriers,  here  assumed  independent  of  E, 
and  the  integral  is  over  the  volume  occupied  by  the  field.  The  prefactor 
B  can,  in  general,  depend  on  E  since  the  electron  energy  distribution 
changes  with  field.  However,  for  simplicity,  this  effect  is  neglected. 


To  perform  the  integral  in  Eq.  (17)  we  must  assume  some  form  for 
the  spatial  distribution  of  E.  First,  we  take  the  radial  dependence  of  the 
field  to  be  gaussian: 


E  =  E  (p,t) 


EQf(t)  e 


(18) 
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Here  f(t)  is  a  normalized  pulse  shape  function  and  0  is  a  parameter  related 
to  the  beam  diameter.  The  bar  over  E  and  EQ  signifies  an  average  of  the 
field  magnitudes  taken  over  one  cycle  of  the  optical  frequency  during 
which  f(t)  is  so  slowly  varying  as  to  be  considered  constant. 


The  second  assumption  relates  to  the  variation  of  E  in  the  depth 
or  z  dimension.  Here  we  assume  that  there  exists  some  characteristic 
dimension  zQ,  such  as  the  focal  depth,  ^he  distance  below  the  surface 
in  which  damage  can  be  detected,  or  a  distance  over  which  n  is 
substantially  greater  than  in  the  bulk  due  to  surface  states.  Then  we 
assume  T?  to  be  constant  for  0  z  £  z^  and  to  be  zero  otherwise. 

Finally,  the  probability  per  unit  time  of  a  lucky  electron  reaching 

e  i  is 

h(t)«  nB  y  exp  [  -K/(EQ{(t)  e~^p  )]  2tt  pdpdz 

»  7rzcnB  ^exp[  -K(l+  /3  p2) /TT^fft)  ]  d(p^)  . 

The  second  step  is  valid  here  if  K fEQ  >  >  1  so  that  Pp2  is  a  small  quantity 
anywhere  that  the  exponential  factor  is  appreciable. 


Performing  the  indicated  integrations,  we  obtain 

Ent(t) 

h(t)  =  7rzcnB  (  )  exp  (-K/E0f(t)  ) 


(20) 


The  interpretation  of  this  result 
effective  field  volume  by 

z  TLf(t) 

Veff  =  *  T  “IT- 


is  aided  by  defining  a  time-dependent 


(21) 


In  other  words,  Eq.  (13)  is  still  valid  for  the  time  dependent  case,  but  with 
N  replaced  by  nV  This  comes  about  because  at  small  fields'.,  o  in  ma¬ 
terials  with  large  K,  the  exponential  factor  varies  so  rapidly  with  E  that 
only  that  portion  of  the  volume  near  the  very  peak  of  the  field  spatial 
distribution  participates  in  breakdown. 
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Now,  h(t)  is  still  not  the  quantity  to  be  compared  to  the  streak 
camera  data,  since  we  have  not  allowed  the  possibility  that  destructive 
damage  occurs  at  some  time  before  t.  Let  us  suppose,  as  discussed 
in  the  previous  section,  that  damage  occurs  with  high  probability  after 
some  electron  has  been  accelerated  to  Further  assume  that  this 
damage  occurs  with  a  very  small  time  delay,  say  less  than  1  nsec. 

This  latter  assumption  is  plausible  from  the  work  of  Seitz,  who  showed 
that  only  some  40  generations  of  an  avalanche  are  needed  to  disrupt 
a  lattice.  In  our  case  the  first  generation  takes  only  M  cycles  of  the 
field,  and  subsequent  ones  may  take  even  less  if,  as  discussed  previously, 
the  resulting  carriers  from  an  ionization  event  have  much  more  than 
thermal  energy.  Therefore,  the  time  delay  between  the  occurrence 
of  the  lucky  electron  and  visible  damage  is  expected  to  be  only  a  few 
hundred  cycles  at  most. 


Consider  an  ensemble  of  samples,  each  subjected  to  the  same  laser 
pulse.  Let  g(t)  dt  be  the  fractional  number  of  samples  in  the  ensemble 
which  damf  ge  between  t  and  t+dt  and  v  (t)  be  the  fractional  number  which 
have  not  damaged  by  time  t.  The  time  origin,  t  =  0,  is  taken  to  be  before 
the  application  of  the  laser  pulse  so  that,  v  (0)  s  1.  g(t)  must  satisfy 

g(t)  -  -  =  v  (t)  h(t)  (22) 


since,  by  our  assumptions,  h(t)  is  now  the  damage  rate.  Therefore 

t 

-I 


v  (t)  =  exp 


h(t ' )  dt' 


L  0 


(23) 


and 


g(t)  =  M*)  exp 


L 

I 


h(t' )  dt' 


(24) 
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The  moot  probable  time  for  breakdown  to  begin  is  the  time,  t^j, 
when  g(t)  is  maximum.  Thus  we  consider 

% 

(^-  <»<V2)  -pf-jMfXt')  =0  . 


(25> 


Since  exp  (-  J  h(t ')  dt' )  *  0,  Eq.  (25)  ho’ds  when 


d(h(tM))  _  ,uu  J 
~m - “  wl  * 


(26) 


Between  t  =  0  and  tEM>  the  time  where  the  field  strength  is  largest 
dh(t)/dt>  0  and  Eq.  (26)  can  be  solved.  For  times  greater  than  tEj^, 
d(h(t))/dt  <  0  and  there  is  no  solution  to  Eq.  (26).  We  conclude  therefore 
that  tM  <;  t^M  which  is  in  agreement  with  the  experimental  results  in 
Figs.  30  and  31.  Note  that  this  result  is  independent  of  the  precise 
form  of  h(t).  The  only  requirement  on  h(t)  is  the  physically  reasonable 
condition  that 


dh 

HE 


>  0  for  all  t  . 


The  formalism  developed  in  this  section  will  be  used  for  detailed 
computer  computations  of  the  temporal  distributions  in  the  rest  of  the 
program. 

C .  Dependence  of  Damage  on  Optical  Frequency 


If  the  laser  wavelength  is  changed,  the  formalism  given  in  the 
previous  sections  would  predict  only  a  change  in  K  (see  Eq.  (15)).  How¬ 
ever,  the  prefactor  A  may  also  have  a  wavelength  dependence  if  the 
initial  carrier  density  n  is  due  to  multiphoton  absorption  of  the  laser 
light  rather  than  thermal  excitation.  Then  A  will  also  depend  on  E, 
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i.e. ,  if  j  photons  are  required  to  excite  carriers  from  the  valence 
band  or  an  impurity  level,  A  will  be  proportional  to  E2K 

This  possibility  will  be  explored  further  when  more  data  on 
damage  probabilities  are  available  using  ruby  laser.  An  attempt 
will  be  made  to  find  a  physically  reasonable  value  of  j  for  the  various 
materials  which  will  be  consistent  with  the  probability  vs  field  data, 
its  wavelength  dependence  (see  Figs.  25  and  °9)  and  the  temporal 
distributions  (see  Figs.  30  and  31). 
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IV .  THEORETICAL  ANALYSIS  OF  SELF-FOCUSING 


The  destructive  effects  of  a  high-power  laser  beam  upon  an  optical 
crystal  will  depend  upon  the  electric  field  intensity  within  the  solid.  How¬ 
ever,  because  of  the  self-focusing  of  the  light  beam,  the  electric  field  within 
the  solid  will  vary  from  point  to  point  in  a  manner  which  must  be  determined 
from  the  intensity  and  the  phase  distribution  of  the  incident  radiation.  To 
provide  a  tool  for  calculating  this  intensity,  we  have  developed  a  computer 
program  which  numerically  calculates  the  beam  dimensions  and  the  electric 
field  in+ensuy  within  the  solid,  given  certain  input  information  which  includes 
the  entrance  face  conditions  for  the  beam,  the  pulse  shape  as  a  function  of 
time,  and  a  formula  relating  the  index-of- refraction  variations  to  the 
variations  in  electric  field  intensity. 

We  have  taken  advantage  of  certain  approximations  which  are  normally 
valid  for  the  laser  beams  of  interest,  in  order  to  simplify  the  mathematics  and 
also  to  reduce  the  amount  of  data  which  is  required  to  specify  the  solution  to 
our  computation.  For  example,  we  make  the  paraxial  approximation,  that  is, 
we  assume  that  the  distance  along  the  beam  in  which  significant  changes  take 
place  in  the  tiv e-averaged  intensity  is  large  compared  with  the  radius  of  the 
beam.  Also  we  assume  axial  symmetry  and  show  that  these  approximations 
are  sufficient  to  lead  to  a  self-similar  solution  near  the  axis.  By  self-similar, 
we  mean  a  solution  in  which  the  power  flux  has  the  form 

I  =  Iq(z)  h2  (?)  ,  (27) 

where  ?  =  r/ a(z),  and  I  and  a  are  independent  of  r  (we  have  suppressed  a  time 
dependence  which  occurs  over  a  scale  comparable  to  the  pulse  duration,  on  the 
order  01  nanoseconds  in  the  experimental  setup). 

Our  mathematical  methods  are  closely  related  to  those  described 
previously  in  Ref  ;-  28-32,  although  these  reference^  generally  refer  to  Gaussian 
distributions  of  intensity,  instead  of  the  more  general  seif-similar  solutions. 
Self-similar  solutions  have  been  used  effectively  in  the  past  for  simplifying  the 
partia1  differential  equations  of  Hydrodynamics;  these  self-similar  solutions 
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might  be  appropriately  referred  to  as  "automodel  flows",  as  in  the  Russian 
hydrodynamic  literature,  since  the  intensity  distribution  across  the  beam  at 
one  position  is  modeled  after  its  own  shape  at  a  different  position.  33-35 

There  will  generally  be  some  response  time  for  the  nonlinear 
response  of  the  dielectric  constant  to  the  applied  field.  For  example,  if  the 
induced  change  of  dielectric  constant  were  due  to  thermal  effects,  there  would 
be  a  time  lag  before  the  temperature  would  build  up  to  its  final  value,  and 
similarly,  after  the  main  pulse  was  turned  off,  a  temporary  trail  of  an  induced 
mhomogeneity  in  the  index  of  refraction  would  remain  and  might  provide  lensing 
effects  during  times  comparable  to  the  response  time.  At  the  oresent  time, 
we  have  incorporated  "memory"  effects  into  the  computer  program  by  assuming 
that  the  index  of  refraction  at  a  point  depends  upon  the  entire  history  of  the 
electric  intensity  at  that  point.  Later  it  would  be  useful  to  include  the  effects 
of  elastic  waves  induced  by  electrostrictive  forces,  which  differ  from  the 
present  model  because  the  acoustic  disturbance  at  a  point  depends  on  the 
history  of  the  electromotive  forces  at  all  points  within  some  neighborhood 
of  the  point  of  interest. 


Dependence  of  the  Index  of  Refraction  Upon  the  Electric  Field 


In  most  of  the  literature  on  self-focusing,  the  index  of  refraction  has 
been  written  in  the  form 


11  =  ^o  +  ^2  U "  ^4  u  ’  (28) 

where  u  =  1/2  t\q  |e2|  with  E  the  peak  electric  field  (E  =  iE  e1Ut)-  Equation  (28) 
is  valid  when  p  changes  slowly  during  a  period  comparable  with  a  period  of  the 
optical  field  but  changes  rapidly  compared  with  the  laser  pulse  duration. 


Equation  (28)  should  be  looked  at  as  the  first  terms  of  a  Taylor' s 
expansion  in  powers  of  the  square  of  the  electric  field.  For  very  strong  fields, 
the  expansion  will  not  be  valid  and  it  is  known  that,  with  a  sufficient  increase  in 
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the  field*  the  mdex  of  refraction  will  approach  a  limiting  value  known  as  the 
saturation  fact,  the  terminEq.  (28)  involving  ^  was  included  by 

Akhmanov  et  al.  as  a  first  approximation  to  the  polarization  saturation, 
ui  general,  the  index  of  refraction  is  some  arbitrary  function  of  E2, 
n  =  f(E  ),  and  our  computer  program  was  written  in  such  a  way  that  any 
desirable  formula  can  be  included  with  little  difficulty.  For  example, 
Wagner  et  al.  0  use  a  relationship  of  the  form 


V’2 

V  °  1  t  Ez/Ez  •  (29' 


where  Eq.  (29)  is  slightly  different  in  form  from  Ref.  16  because  the  nonlinear 

effects  in  that  reference  were  expressed  in  terms  of  dielectric  susceptibilities 
instead  of  the  indices  of  refraction. 


The  effect  of  relaxation  processes  have  been  included  in  our  computer 
program  by  generalizing  the  definition  of  u  in  Eq.  (28)  to 


u  = 


2 


(30) 


where  tq  is  a  relaxation  time.  The  change  in  E  in  Eq.  (30)  will  occur  over 
periods  comparable  to  the  pulse  duration  which  in  our  experiment  will  be  on 
the  order  of  nanoseconds.  If  the  relaxation  time  is  very  short  compared  to  the 
pulse  duration,  Eq.  (30)  reduces  to  the  instantaneous  form 


u 


1 

2 


"o  Ie2<‘> 


(31) 


which  was  given  previously.  On  Ihe  other  hand,  when  the  relaxation  time  is  long, 
then  Eq.  (30)  approaches  the  form 


where  the  pulse  has  been  taken  to  vanish  for  t  <  0. 
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(32) 


More  realistically,  when  the  relationship  between  the  index  of  refraction 
and  the  field  is  a  linear  functional  of  E  ,  we  would  have 

t 

u(r,t)  =  \  j*  dT  Jdr'  G(r",t")  |e2(t;P)  |  ,  (33) 

“OC 

where  t"  =  t  -  Tand  P1  =  r  -  P .  Equation  (33)  is  the  sort  of  relationship  we 
might  expect  if  heat  or  momentum  were  released  at  a  position  P  at  time  Tand 
G  plays  a  role  similar  to  a  Green1  s  function  in  describing  the  propagation  of  the 
effect  to  the  position  r  at  time  t.  We  now  make  the  somewhat  drastic  approxi — tion 
that 


G(r"(t")  =  6(?V  g  (t") 

which  ignores  the  transfer  of  fluctuations  in  the  index  of  refraction  from  point  to 
point  in  the  medium  as  would  be  expected  to  occur  due  to  thermal  conduction  and 
convection  when  the  fluctuations  are  caused  by  thermal  effects  or  due  to  elastic 
waves  when  the  fluctuations  are  caused  by  electrostriction.  Then  Eq.  (33)  becomes 

t 

u(r,t)  =  j  tio  y  '|E2(t,P)  |  g(t")  dr  .  (34) 

-a 


Equation  (30)  results  from  Eq.  (34)  when  we  take  g(t")  =  exp  (-t"/ T0  )/t 
which,  with  an  appropriate  choice  of  the  relaxation  time  t  ,  should  be  an  appropriate 
approximation  since  it  is  physically  apparent  that  g(t")  must  be  a  monotonically 
decreasing  function  of  t". 

E  •  Electrostrictive  Effects 

The  characteristic  time  scale  for  elastic  v. aves  excited  by  a  laser  beam 
is  on  the  order  of  the  beam  radius  divided  by  the  longitudinal  sound  velocity. 
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Although  the  precise  value  depends  upon  the  beam  radius,  this  suggests  a 

“8 

characteristic  time  on  the  order  of  10  seconds.  Since  this  is  comparable 

to  the  pulse  duration,  and  also  because  the  exit-face  surface  damage  might 

involve  elastic  waves  striking  a  semimolten  surface,  one  might  wish  to  include 

these  effects  in  the  final  version  of  the  computer  program  as  a  separate  term 

independent  of  the  relaxation  form  described  in  the  previous  section,  using  an 

32 

approach  similar  to  that  of  Kerr. 


The  electrostrictive  force  is  a  volume  force  proportional  to  the 
gradient  of  the  intensity,  i.e.  , 


f  =  aVl 

12  _ » 

where  I  =  ^-E  .  If  we  let  u  be  the  displacement,  the  equation  of  motion  is 

pu  =  F  +  f  ,  (35) 


where 


9  a 


Fi  = 


'ik 


For  a  solid  satisfying  Hooke's  law,  the  stress  is  a  linear  function  of  the  strain 

3  6 

and  when  the  body  is  also  isotropic,  we  have 

u  =  c^  grad  div  u  -  c^  curl  curl  u 

+  grad  i  ,  (36) 

P 

2  2 

where  c^  =  E(1  -  a)/p(Ho),  =  E /  2p  (l+o) ,  E  is  Young' s  modulus ,  and  a  is 

Poisson1  s  ratio. 


Since  an  arbitrary  vector  can  be  written  in  the  form 
u  =  -  grad  curl  A  ,  (37) 
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Eq.  (36)  becomes 


grad  V  <£  +  —  ij  =  curl  £  A  +  c2  curl  curl  A  J  .  (38) 

In  an  infinite  medium,  it  is  apparent  from  Eq.  (38)  that  only  the 
longitudinal  wave  will  be  excited  by  electrostrictive  forces  and  that  it  will 
satisfy 

2  2  a 

l  =  c£  V  ^"p“  1  ‘  (39) 

At  first,  we  will  neglect  the  vector  wave  A*  although  we  note  that  trans¬ 
verse  waves  described  by  the  vector  A*  will  be  excited  by  boundary  conditions  at 
the  surface.  This  effect  may  be  important  for  calculations  of  the  acoustic  intensity 
at  the  exit  face. 


Now  let  p'  =  p  -  pQ  be  the  increase  in  density.  From  the  continuity 
equation  for  small  amplitude  waves,  we  find 


so  the  wave  equation  for  p'  becomes 


3  t 


a  V2  I 


Equation  (40)  is  the  wave  equation  with  a  source  term  depending  upon  the 
optical  intensity  I.  The  transverse  Laplacian  terms  in  a  paraxial  laser  beam  will 
be  much  larger  than  the  longitudinal  terms.  Then,  Eq.  (40)  can  be  solved  as  a 
two-dimensional  problem  at  each  point  z,  and  for  axial  symmetry  has  the  form 


9  V  -  2  I  JL  r-  <L£  -  £  9  „  9  i 

7 . 2  X  r  3  r  *  9  r  r3rr3r 


(41) 
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Equation  (41)  can  be  solved  numerically  by  making  use  of  the  usual 
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Green1  s  functions  by  methods  similar  to  those  of  Kerr  and  must  be  solved 

simultaneously  with  the  equation  for  the  intensity  I.  When  the  laser  pulse  is 

turned  on,  Eq.  (41)  describes  a  sound  wave  which  initially  implodes  upon  the 

axis,  giving  a  high  density  on  the  axis,  thus  forming  a  self-focusing  acoustic 

lens.  As  time  goes  on  the  density  on  the  axis  will  temporarily  decrease  as  the 

acoustic  wave  is  reflected  from  the  axis,  in  a  manner  dependent  upon  the  laser 

pulse  shape,  but  it  will  essentially  reach  the  steady  state  condition  in  which  the 

left-hand  side  of  Eq.  (41)  vanishes.  Then  p'  takes  the  steady-state  form 

P*  =  —g—  I  .  (42) 


with  an  index  of  refraction 


(43) 


Thus  for  times  large  compared  with  the  relaxation  time,  the  steady- 
state  form  of  Eq.  (43)  is  the  same  form  as  Eq.  (28)  (with  =  0). 

C  The  Paraxial  Beam-Tracing  Equation 

We  consider  a  laser  beam  propagating  in  the  z -direction  and  symmetrical 
about  the  z-axis.  We  assume  that  the  electric  field  magnitude  can  be  written  in 
the  similarity  form 


E(r ,  z)  =  EQh(£)a(o)/ a(z) 


(44) 


where  f  =  r/  a,  Eq  is  independent  of  r  and  z  and  a  is  a  function  of  z  alone.  The 
radial  distribution,  h(£),  is  determined  by  the  boundary  conditions  at  the  entrance 
plane,  and  may,  for  example,  have  the  Gaussian  form 


h(£) 


e-<€2/2> 


(45) 
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Here  a  may  be  considered  as  the  beam  radius.  We  have  suppressed  a  slow 
time  dependence  in  Eq  and  a,  which  change  over  a  time  comparable  with 
the  pulse  duration. 

Although  we  have  assumed  the  self-similar  form  given  in  Eq.  (43), 

we  will  show  that,  near  the  axis  of  the  beam,  this  form  is  consistent  with 

the  governing  equations.  These  governing  equations  can  equally  well  be 

2  9 

obtained  from  the  ray -tracing  equation  of  geometrical  optics  or  from  the 

28  30 

eikonal  approximation  to  the  wave  equation.  '  In  the  former  approach, 
the  diffraction  term  must  be  added  in  an  ad  hoc  manner  to  give  agreement 
with  the  linear  wave  equation. 


Taking  the  wave  equation  approach,  with  vanishing  conductivity, 
the  wave  equation  is 


curl  curl  E  = 


a  2d 


at 


(46) 


/ 


where  we  have  assumed  the  magnetic  permeability  p  =  1.  The  electric  dis- 
— »  2  — # 

placement  D=  t)  E  where  ti  is  the  index  of  refraction  and  is  inhomogeneous 

2 

since  it  depends  upon  E  .  In  general  we  assume  that  r\  depends  upon  the 
2 

average  of  E  over  many  optical  cycles  so  that  to  sufficient  accuracy 


a  2d 
a  t2 


~  2  a2E 

=  n  — t- 


at 


and  Eq.  (46)  becomes 


V 


_  2 
grad  div  E  = 

c 


32E 


at 


(47) 


Now 


o  »  ...  —4 

div  D  =  0  =  t\  div  E  +  2  r^E  *  V  r|  , 
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so  div  E  =  -  E  *  V  (  In  (rf)).  Thus  the  div  if  term  does  not  generally  vanish; 

32 

however,  it  is  generally  small  except  near  the  focal  point  and  we  will  neglect 

30 

the  term  henceforth.  To  proceed  we  assume  that 
£.•§!<!>  Eo(t)h(S)e-iks<r’2)eI(ut-kz) 

=  e  4.E  «)  ei(ut'kz)  ,  (48) 

T  o 

2  2 

where  we  specify  EQ(t).  h{§)  and  s(r,0)  =  r  /  2R  at  the  entrance  face  where 
the  boundary  condition  on  the  eikonal  s(r,0)  has  been  chosen  so  that  the  light 
will  focus  at  the  point  z  =  R.  The  propagation  constant  k  =  c* 


Neglecting  div  E,  Eq.  (47)  becomes 


i  JL  r  <14. 

r  9  r  9r 


2  ,  i/ 

"  Vt 

c 


ip  =  0 


(49) 


g  u  i 

For  paraxial  beams  the  term  — {jp-  can  be  assumed  small  in  comparison  with 
the  transverse  Laplacian.  Th£nzlet  ^  ^  e’iks,  multiply  Eq.  (49)  by  eiks 

and  equating  imaginary  parts ,  we  obtain 


1_ 

r 


a 

ar 


r  & 


a a 
ar 


(50) 


Integrating  Eq.  (50)  over  a  volume  bounded  by  r  =  R  and  the 
pianes  z  -  z^,  and  z  ~  Zg,  when  R^  *oo,  shows  that  the  power 

R 

2ir  J  ^o  r  dr 
0 
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is  a  constant  independent  of  z  which  is  necessary  in  a  nonabsorbing  medium. 
This  is  consistent  with  the  self  -similar  form  since 


r  dr 


R 

j*  h2(?)  r  dr 
0 

R/a 

J  h2(?)  ?  d? 
0 


constant. 


In  general,  Eqs.  (50)  and  (51/  are  in  a  convenient  form  for  numerical 
solution  independent  of  the  similarity  assumption,  since  these  equations  could 
be  numerically  integrated  by  a  marching  method  starting  from  the  boundary 
conditions  at  z  =  0.  For  a  similarity  solution  we  expand  s,  and  t|  in  Taylor's 
series  in  f ,  which  have  the  form 


s(r ,  z)  =  sq(z)  +  i  s2(z)  ?2 


K0U)  +  ‘  K2(z)  52 


(53a) 

(53b) 
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so  that  substituting  \Eq.  (53a)  into  Eq.  (50)  yields 


da 

a  Hz  =  s2 


\ 


(54) 


while  substituting  Eq.  (ii3b)  into  Eq.  (52)  yields 
9s^  o  Id  s  \  2s-  ,  s?  \ 

2  sr  ♦ «  [t#  '*■  nr  £  *  4  -  Ko<z>  - 
\  a  / 

1  1  9  .  9h 

h5  8T  ?  9? 


K2(z)f-‘ 
- 2 — 


k  a 


(55) 


If  h  is  the  Gaussian  e‘?  /2,  the  right  hand  side  of  Eq.  (55)  is  (?2  -  2)/k2a2. 
For  other  functions  we  can  expand  d=l/h£9/9£(£  9h/3$)ina  Taylor 1  s 
series  so  that  the  right  hand  side  equals  d(0)  +  d"  £2/2.  In  either  case  we 
obtain  an  equation  for  s2  by  equating  the  coefficients  of  £2.  For  the  Gaussian 
we  obtain 


d 

~ST 


S2  + 


a  (z) 


K2(Z)  j 

?? 


(56) 


while  an  equation  for  sq  can  be  obtained  by  equating  the  constant  term.  (The 
axial  eikonal  sq  is  not  needed  to  determine  the  intensity  but  it  could  be  useful 
in  determining  the  spectral  output  of  a  nonlinear  beam. )  Substituting  Eq.  (56) 
into  Eq.  (54)  gives  the  following  second-order  equation  for  a, 


Equation  (57)  can  be  considered  the  equation  for  the  radius  of  the 

laser  beam.  For  a  Gaussian  beam  the  power  in  any  cross-section  falls  off 
2  2 

as  exp  (-  r  /a  ).  The  initial  conditions  require  the  beam  radius  a  and  the 

o 

first  derivative  da/dz  =  -aQ/R.  The  total  power  delivered  to  the  solid  is 
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00 


c  a?  E?(t)  j*  r\(V  ?h2(C)  d? 


o  o 


0 


c  2  -,2  ,,v 

tjt  ft  (t)  T) 

o  O  O  O 


(58) 


which  determines  EQ(t),  given  aQ  and  the  instantaneous  total  power  P» 


Whei.  q  is  given  s  a  function  p  =  f(u),  Eq.  (57)  can  be  written  as 

2 


A..  Lu)|£s^i 
a?  I  \  77  77 


(59) 


1  2 

and  when  u  has  the  instantaneous  form  u  =  ■=■  p  E  ,  Eq.  (59)  becomes 

«  o 


j2 

d  a 

JT 

dz 


E2a2 
o  o 


n  a 
'o 


h(0)  h"(0)  f(u) 


^  U  /r- 


n 

0  k  a 


(60) 


where  u  at  r  =  0  has  the  value 

i  '>oEo  (t)  h2,0)  • 

In  the  more  general  case,  the  computer  program  solves  Eq.  (59)  where  u 
is  evaluated  by  knowing  the  value  of  E  at  a  given  point  over  its  entire  past 
history,  and  applying  Eq.  (30). 


D.  Numerical  Methods 


In  the  case  where  the  beam-tracing  equation  is  used  with  an 

30 

instantaneous  response  mechanism,  Wagner  et  al.  have  noted  that  the 
equation  has  the  same  form  as  the  equations  of  motion  for  a  single  particle 
in  a  potential  well,  which  allows  one  to  make  use  of  analytical  solutions  fcr 
certain  relation?  hips  between  the  index  of  refraction  and  the  fieldstrength. 
The  approach  used  in  the  computer  program  is  to  numerically  integrate  the 
equation,  a  procedure  which  allows  us  to  be  flexible  in  choosing  different 
functions  r\  =  f(u)  and  also  to  include  response  time  effects  which  are 
incorporated  in  Eq.  (30). 
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In  the  numerical  computation  the  function  h(f )  has  been  taken  to  be 
a  Gaussian.  The  values  of  u  on  the  axis  at  time  t  are  stored  at  up  to  400 
positions  along  the  center  of  the  beam.  The  increments  in  u^  from  t  to 
t  +  A  t  are  then  calculated  by  a  method  which  amounts  to  a  finite  difference 
expansion  based  on  Eq.  (30).  However,  this  updating  requires  the  values 
of  the  radius  a  at  t  and  t  +  A  t,  that  are  found  simultaneously  by  integrating 
the  beam -tracing  equation  (Eq.  (40))  from  the  entrance  face,  this  being 
accomplished  by  a  fourth-order  Runge-Kutta  integration.  During  this  integra¬ 
tion,  a  test  for  accuracy  is  made  and  the  step  size  in  the  integration  is 
automatically  divided  by  two  so  many  times  as  is  necessary  until  the  criterion 
is  met.  An  alternative  program  was  written  to  use  only  the  instantaneous 
version  of  Eq.  (30)  which  gave  a  convenient  method  of  testing  the  accuracy  of 
the  s  econd  program  by  letting  the  relaxation  time  in  that  program  to  be  taken 
extremely  small  compared  with  the  pulse  duration.  The  instantaneous  program 
and  the  program  including  relaxation  are  described  in  the  next  two  sections. 


E*  Beam-Tracing  Equation  with  an  Instantaneous  Response  Mechanism 

In  order  t  o  b  e  able  to  predict  the  field  intensities  along  the  track, 
two  programs  were  written.  The  first  program  treats  the  case  where  the  index 
of  refraction  along  the  track  adjusts  instantaneously  to  the  applied  field.  The 
second  program  which  will  be  described  in  Sec.  IV-F  treats  the  case  where  there 
is  a  time  lag  between  the  application  of  the  optical  field  and  the  adjustment  of  the 
index  of  refraction  along  the  trail.  The  first  program  is  quicker  to  run  and  a 
can  be  operated  with  a  larger  step  increment  in  the  z -direction. 


When  relaxation  effects  a  re  ignored  the  mathematical  problem  can  be 
expressed  in  the  form  of  Eq.  (60)  which  can  be  rewritten  as 


where 


U 


h"(0)  (  f2 

shut  pr 


+ 


1 


Tnn 

2k  a 


(0. 


\ 

/ 


(62) 


r=  0 
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where  (for  a  Gaussian  beam  satisfying  Eq.  (45)  so  that  h"(0)  =  -  h(0)) 

Eqs.  (61)  and  (62)  are  equivalent  to  Eq3.  (2.16)  and  (2.17)  of  Wagner  et  al.  , 30 
except  for  an  unimportant  constant  in  the  definition  of  U  and  a  slight  change  of 
notation.  (We  use  f/  t)q  for  the  ratio  of  the  nonlinear  diffraction  index  to  the 
undisturbed  index  while  they  iefine  e  ^  +  X  as  the  susceptibility  with  constant 
magnetic  permeability  so  that  1  +  x  /  c  =  f 2/r)^. ) 

If  we  let  p  =  da/ dz,  Eq.  (61)  is  easily  integrated  to  give 

1  2  TT  „ 

2  p  U  =  E  (63) 

where  E,  "the  energy",  is  a  constant  as  we  move  along  the  beam.  It  can  now 
be  easily  recognized  that  Eq.  (61)  is  equivalent  to  a  particle  moving  in  a  one¬ 
dimensional  c  onservative  potential  U,  as  noted  by  Wagner  et  al.  ,  30  or  to  the 

O  O 

two-body  central-force  problem  where  U  is  the  fictitious  potential,  provided 
that  we  interpret  z,  the  distance  along  the  beam  as  time,  and  a,  the  radius  of 
the  beam,  as  position,  with  p  =  da/dz  playing  the  role  of  momentum.  The 
identification  of  our  problem  with  the  standard  problem  of  classical  mechanics 
allows  us  to  use  our  familiarity  with  these  solutions  in  the  beam-tracing  problem. 
Wagner  et  al.  30  give  a  few  illustrations  of  this.  Another  illustration  is  the  "giant 
planet"  instability;  that  is,  it  is  known  that  a  relatively  small  perturbation  caused 
when  a  comet  passes  near  a  large  planet  may  change  the  orbit  in  such  a  way  that 
a  previously  captive  comet  may  be  thrown  into  an  escape  orbit.  A  similar  effect 
happened  in  the  original  version  of  beam -tracing  program  in  which  we  found  that 
small  numerical  errors  would  cause  the  escape  of  a  beam  which  had  previously 
been  trapped  and  had  described  several  focal  points;  the  inverse  effect  was 
probably  possible  in  which  a  numerical  error  could  cause  a  beam  to  be  trapped 
even  though  its  initial  conditions  should  not  have  allowed  this.  This  nurrerical 
problem  was  solved  in  the  final  version  of  the  program  by  requiring  the  derivative 
da/ dz  to  satisfy  the  constraint  represented  by  Eq.  (63). 

Another  illustration  concerns  an  expression  in  common  usage  for  the 
nonlinear  index  of  refraction 


2 

n 


(64) 
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which  leads  to  the  potential 


u  =  -i  (1  +  A 

*  n 

'o 


E2  a2 

-V->  + 


2??" 


or  equivalently  by  adding  a  constant  to  U,  we  replace  U  by 


U' 


2k  a 


Y~T  (1  k2  E2  a2  ) 

4-*  n  O  o 


(65) 


(Note  that  in  Eq.  (64)  we  have  departed  from  the  general  form  given  by 
Eq.  (28).  Unfortunately,  there  is  a  great  deal  of  confusion  about  these  expansions. 
For  example,  Akhmanov  et  al.  state  that  they  are  expanding  the  index  of 
refraction  as  rj  =  +  r?2  |  A  | 2  +  r?4  |  A  |  4  but  when  they  apply  ^  ^  their 

derivation  of  their  Eq.  (8)  and  thereafter,  they  have  actually  used  the  expansion 


(66) 


This  means  that  the  usual  form  keeping  only  the  quadratic  terms  in 
Eq.  (66)  involves  a  fourth  order  term  in  the  index  of  refraction.  All  this  means 
is  that  Akhmanov  et  al.  have  really  expanded  in  terms  of  susceptibility  like 
Wagner  et  al.  and  there  is  really  no  physics  lost,  but  the  inconsistency  can  be 
confusing.  In  any  case,  our  use  of  Eq.  (64)  allows  our  results  to  be  directly 
comparable  to  Eq.  (18)  of  Akhmanov  and,  by  replacing  n«/rj  by  /e  ,  to 
Eq.  (3. 2)  of  Wagner  et  al. )  o  ^  l 

Now,  it  is  immediately  apparent  from  the  analog  with  mechanics  that, 

if  the  initial  conditions  at  z  =  0  are  a  -  a0  and  da/dz  .  a'  ,  the  beam  will  be 
trapped  only  if 
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in  which  case  there  are  an  infinite  number  of  focal  points.  On  the  other  hand, 

in  the  opposite  case  there  will  either  be  no  focal  points  or,  at  most,  one  (if 

28 

a^/  aQ  <  0).  Now  Akhmanov  et  al.  state  conditions  in  which  they  claim  vhere 
are  exactly  two  focii.  Actually  they  have  introduced  a  spurious  solution  in  going 
from  their  Eq.  (19)  toEq.  (20).  Their  Eq.  (19)  has  the  form 

<§>2  ■'£  +  «,.  (68) 

where  c ^  and  are  constants  defined  in  their  text. 

An  integration  of  Eq.  (68)  gives 


z  h  -r2  Vci +  c2  3  ^  Vci +  c2 f2  •  (6S 

with  the  boundary  condition  f  =  1  at  z  =  0.  Now,  Akhmanov  et  al.  give  as  the 
integral  of  Eq.  (68)  the  result  of  squaring  Eq.  (69),  that  is 


c2z 


2  +  2  (Cj  +  c//2 


z  +  1 


(70) 


which  is  also  equivalent  to  Eq.  (3.4)  of  Wagner  et  al. 

At  a  focal  point  f  =  0,  so  Akhmanov  et  al.  assumes  that  the  values  of 

z  at  which  focal  points  occur  is  obtained  by  setting  the  quadratic  on  the  right- 

hand  side  of  Eq.  (70)  equal  to  zero,  which,  for  a  positive  discriminant,  leads  to 

two  solutions.  However,  one  of  these  solutions  is  spurious  and  only  the  solution 

corresponding  to  setting  f  =  0  in  Eq.  (69)  can  be  real.  Thus,  the  results  plotted 

in  Fig.  1  of  reference  28  are  incorrect.  This  error  would  be  difficult  to  spot  if 

we  did  noi  hove  guidance  from  the  analogy  with  classical  mechanics.  (Since  we 

28 

have  gone  to  some  trouble  pointing  out  small  mistakes  in  Akhmanov  et  al.  ,  we 
should  state  that,  taken  as  a  whole,  their  paper  is  not  only  useful  but  also  one  of 
the  important  basic  references  in  the  theory  of  self-focusing.) 

A  listing  of  the  program  which  was  written  to  integrate  the  beam  tracing 
program  is  given  in  Fig.  32  a-d.  The  shape  of  the  pulse  into  the  medium  can  be 
chosen  arbitrarily  by  adjusting  constants  in  the  DATA  statement  in  subroutine 
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PBN-506 


C  DTrTIME  SA  VP  LE  .  N  YE  T=  NO  OF  FRAMES 

C  THIS  PROGRAM  CALCULATES  THE  SELF  FOCUSING  OF  A  LASER  BEAM 
COMMON  HEX.CAP.H.DEE 

DIMENSION  A!  401)  .P  1401  >.Z(  401)  .POW  (12) 

NAMELIST /PARAM/  RELAX. E2.E4 
READ  (5  .PARAM  ) 

WRITE!  6.  PARAM) 

NAME  LI  ST  /CONTR  0/  D  T  .  NY  ET  .L  EG  UP  .N  A.  NB  .LETS.  NET.  INTO 
READ  (  5.  CONTRO  ) 

WRITE!  6. CONTRO) 

60  0  RE  AD  (5  .100)  G.  WL.R.AL.  AO 

100  FORMAT  (5E10.4) 

C  POWER  IN  MEGAWATTS.  COEFFICIENTS  IN  CM /STV  T.  WVLE  NG  TH  INCM.RAD  OF  CUR 
C  IN  CM,  LENGTH  IN  CM. NO  OF  POINTS  .R  AD  IUS  OF  BEAM  IN  CM 
PCRITr  .0  0  30*  WL  *WL/64./  3.  14  159/  3.  14  159/E2  *2  . 

HrAL/F  LOAT  (NET  ) 

NET  1  rN ET  •*  1 
GIM-  G 
T-0  •  0 

IT  IM=N  YE  T*DT 
WRITE!  6.  903)  ITIM 

90  3  FORMAT  (1  HI. 2SX»*  TIME  SEQUENCES  FOR  A *. 13 »  *  NANOSECOND  PULSE  V  > 
DO  5  I-1.NET1 

2!  I)  =H*FLOAT  (1-1  )  ♦0.1E-06 

A(I)rO.O 
5  P(I)rQ.O 

DO  190  L  OC  K-  1 .  NY  ET 
T-T*DT 

CALL  POWER(T.PAT) 

G=GIM*PAT 
A  (  1 )  rA  0 
P(  1)  =-  AO/R 

CAP=  WL**2/4./3.  14159/  3.  14  159-  8.*G/0.  00  3  *E  2 
HE  Xr  64.  *E4*  (G/0  .003  )**2 


Fig.  32a  Computer  Program  to  Integrate  the  Beam  Tracing 
Problem. 


93 


PBN-507 


OEE=  (A  O/RJ  **  2+  <C  AP  ♦  HEX/Z. /AO/AO  )/AO/A 
nIXl 35 

00  15  J=!»NET 
IF  f  M IX  .N  E.  35  )  GO  TO  83  2 
IF  (A  !J  ).  LT  .0.0  )  HIX  =  J 
IFIIP{  J)  /R)  .GT.  0.0)  HIXrj 
83  2  CONTINUE 
ANrAtJ  ) 

PNrP(j> 

CALL  RATEtAN.PN) 

A I  1 )  ZAN 
15  P!  J*l)  r  pn 

IF!MOD!LOCK*LETS).NE.O)  GO  TO  19  0 

PCWIT#G,£2»£4#WL»R*AL  »A  0*  NE 
210  FORMAT  !  35X.  16HCRITICAL  POWER  r,E12.4*2H) 

2  E12.4.2HMW/ 

3  1  OX  1 1  OX  1 5HN2/  NO  *10X  .5HN4/N0  » 5  X » 

J  1QHWA VELENGTH* 

5  9X.6HEXTENT.  6X  *9HR AD IUS  A0.11X. 

Glz_0  «ihpesh.  /it 

00  413  MrNA.NB  . LEGUP 
LZ-LZ>  1 

*»13  POW!  LZ  l-G/3.  14  159/ A  (M)  ••  2 

WRITE! 6. 200) T, !Z !N) ,NrNA *NB  .LEGUP).  (A  ( 
!P  f  N )  » N-NA  »NB  .LEGUP),  fP 
200  FORMAT  !26X  •  ’BEAM  GEOMETRY  WHEN  Tr  *  .F  6 

2  J «.  .  ,m  m. 


MW/37X  »  14HACTUAL  POWER  r, 

GX,  9H CURVATURE, 
ox  .6  El  5.  5.115/) 


uihi.nzna.NB  *LEGUP).  f  a 
IP  f  N  )  »  N-NA  *NB  tLEGUP)*  f 
•BEAM  GEOMETRY  WHEN  Tr  *  ,F 

/4  X8  HP 


*«■  *  .  ...  -  .  ..  »  “*  *o  nr 

3S*  1 1  El  0.  4/7X  .5HSL0PE  *1  IE  10  .4  /4  X*  8HP0  YN TI 
LZ  =  MIX-€ 

00  199  Nr l.U 
LZrLZ^ 1 

19  9  ROW  (  M )  r  G/3.141  59/A  (LZ)  **2 
LZr  MIX  -5 


(N)*NrNA*NB  tLEGUP). 

POW(N)  »Nr  1  ,11) 

G.  2,  'NANOSECONDS  • 

OSIT  ION*  11  FI  0.  7/6X  .6  HR  AOIU 
NG  *1  IE  10  .4  /4X*  4HFLUX  ) 


LZUMIX*5 
WRITE!  6,  205) 


WRITE!  6*  205)  !Z  !N )  ,  NrLZ  »L  Z1  )*  (A 

2  (P  (N ) » NrLZ  *L  Z1  )*  i 

20  5  FORMAT  !2GX.  ’GEOMETRY  NEAR  THE  BEAM  VA 

m 

°*  <I/7X  *5HSL0PE  *1  IE  10  .4  /4  X*  8HP0  YNTI 
190  CONTINUE 

INTO  r I  NT  0-  1 

IF  !  I  NT  0.  GT  .0  )  GO  TO  GOO 

CALL  EXIT 

ENO 


!N  )*  Nr  LZ  »LZ1  ), 

POW!  N)  *Nr  1  ,H) 

1ST* 

OS  IT  ION.  11FI  0.  7/GX  .6HRA0IU 
NG  *1  IE  10  .4 /4X*  4HFLUX  ////) 


Fig.  32b  Computer  Program  to  Integrate  the  Be*m 
Tracing  Problem. 
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SUBROUTINE  R AT E(  AN  »PN) 

COMMON  HEX.CAP.H.DEE 
F(  Y)r(  CAP*HEX/Y**2)/Y**3 
LAMBrJ 
LETGOr  1 
NOB-1 
OUAO=H 
11  5  CONTINUE 
T 1 1  r  H*  PN 
T1  2-  F|  AN  )  *  H 
T21  =  H*  (PN*0.  5*  T1  21 
BY r  AN*  0.  5*Tl  1 
T22  =  FIBY1  *H 

IFH  II  .-T12/T22)  **2)  -L  T.  .Oil  GO  TO  15 
IF  II  T22/PN)**2.LT.  0.1E-07)  GO  TO  15 
IF  (NOB  .EQ.  -1  .AND  .LAMB. EO  .1  ) 

2  VR  ITE(  6.  139)  LET  GO  •  AN 

IF  I  NOB  .EO.  -1  )  LAHB=-1 

13  9  FORMAT  (IX.  30H  SLIPPAGE  ,  1  3HLETGO  EQUALS  . 

2  E  1 5  .4  ) 

IF  (NOB  .EO.  -I  )  GO  TO  IS 
HrH/2. 

LETGOr  2*  LE TGO 
GO  TO  115 
l  5  CONTINUE 

NOBr-1 

T31  =  H*  (PN*2.  *T 22  -T  12  ) 

BY  =  AN*2.*T2I-T1I 
T32=F(BYi  *H 

ANrAN*  <T11«4.*T21*  T31I/6. 

PNrPN*  (T12+4.  *T22*  T32)/6. 

LETGOr  LE  TGO- 1 
IF  (LET  GO  .GT.  0)  GO  TO  115 
P2rDEE  -(  CAP* HE  X/2.  /AN/ AN  )/ AN/AN 
IF  (P2.  LT  .0.0)  GO  TO  53 
P2rSQRT(  P2) 

PNrP2*ABS(PN)/PN 
5  3  Hr OU AO 
RETURN 
END 


Fig.  32c  Computer  Program  to  Integrate  the  Beam 
Tracing  Problem. 
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20 


300 
40  0 


•  50./.  (FIT  INI.  Nr  1.51  /0..  .7  ..2. 


SUBROUTINE  POWER  (  T.  p*r  } 

DIMENSION  T»U(  IOI.FITI  10  I 

'%Jn  ^1^U,N,  ,Nrl  I/O.  *10.  .20.  .3  0. 

«  .0.0/  •  N/5/ 

-1 

Nr  N*  1 

IFCN.GT.  M)  GO  TO  300 
IF  (TAU  CN  I.  LT.T  »  GO  TO  20 

GO  TO  400*  FIT,N*  ”FIT  <N_l  **  **T-TAU(N-»  1/  (TAU  (N  )~TAU(  N- 1)) 

PATrO.O 

RETURN 

END 


Fig.  32d 


TrTcSg^roM^  t0  Int6grate  the  Beam 
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POWER.  The  peak  power  in  megawatts  is  entered  by  the  READ  statement 
in  the  main  program,  as  is  optical  wavelength  in  cm  inside  the  material, 
the  geometrical  radius  of  curvature  of  the  beam  at  entry,  the  depth  of  the 
sample  to  which  the  calculations  are  to  be  carried  out,  and  the  effective 
radius  of  the  Gaussian  beam.  The  NAMELIST  PARAM  inputs  the  relaxation 
time  in  nanoseconds  and  the  coefficients  for  the  expansion  from  Eq.  (28)  in 
cm/statvolt.  The  NAMELIST  CONTRO  introduces  the  time  increment 
DT  in  nanoseconds,  NET, the  total  number  of  mesh  points,  and  some  additional 
control  integers.  A  typical  example  of  the  input  data  cards  is  shown  in  Fig.  33. 
The  output  which  gives  the  critical  power  for  self-focusing  and  the  actual  power 
in  the  beam  at  a  specified  time,  is  shown  in  Fig.  34.  The  output  gives  the 
radius,  slope  and  Poynting  flux  for  the  beam  at  fairly  widely  separated  positions 
along  the  beam  and  also  at  closely  spaced  intervals  near  the  beam  waist  (or  at 
the  focus  if  the  beam  radius  actually  reaches  zero). 

F*  Numerical  Methods  for  Integrating  the  Beam-Tracing  Equation 
with  a  Finite  Response  Time  Mechanism 

In  general,  there  will  be  a  time  lag  after  the  application  of  the  optical 
field  before  the  index  of  refraction  adjusts  to  its  final  value.  Since  the  perturba¬ 
tion  of  the  index  of  refraction  remains  along  the  path  of  the  beam  for  some  time 
after  the  optical  field  is  reduced,  this  effect  is  also  known  as  the  residual  trail 
effect.  In  our  formulation,  the  temporal  delay  is  given  by  Eq.  (30).  In  the  case 
where  the  index  of  refraction  is  dependent  only  on  the  quadratic  term  in  the  electric 
field,  Eq.  (30)  is,  in  fact,  equivalent  to  an  integration  of  the  differential  relaxation 
equation  given  by  Fleck  and  Kelly39  and  by  Shen  and  Loy.  40  The  expression  for  u 


u(r,  z,t) 


n  a2 
o  o 

~TF~ 


r  Eo(T) 

0  a^(z,t) 


\ 

(71) 


where  we  take  h  to  be  the  Gaussian  distribution  h2  =  exp  (-r2/a2).  The  beam 
tracing  equation  is 
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ax  QT 
1  PA  RA  H 

R  EL  AX  =0  .  I E-  04  . 

E2=0.2E-12#  E4=0.32E-23. 

1END 
1  CO  NT  RO 

0T=2.0.  NYET=13. 

L  EG  UP  -  4  0  »  N  A-  1  •  NB-40I. 

LETS=1.  NE  T"  400# 

INTO:  2. 


SEND 

2  n 

♦  01 

•1060-03 

•25  ♦OO  .5 

♦  00 

.5  00  0-  02 

.211 

♦  01 

•1060-03 

.25  ♦OO  .5 

♦  on 

•5  00  0-  02 

3F  IN 

Fig.  33  An  Example  of  the  Input  Data  to  the  Beam 

Tracing  Program. 
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In  Eqs .  (73)  and  (74)  P(r)  is  the  instantaneous  power,  which  should, 
strictly  speaking  be  replaced  by  P(r  -  z/c).  The  exact  form  is  not  necessary 
in  our  experimental  setup  where  the  thickness  of  the  crystal  L  always  satisfies 

L  «  c  P(t)/P'(t)  ,  (75) 

^Ut  f4°0r  crystals  there  are  some  interesting  effects  as  shown  by  Shen  and 
Loy  in  which  the  focal  spot  can  move  with  a  (virtual)  velocity  which  exceeds 
the  speed  of  light. 


The  mathematical  problem  requires  us  to  solve  the  two  integrals  from 
Eq.  (/3)  and  (74)  simultaneously  with  the  differential  equation,  Eq.  (72),  with  an 
arbitrary  pulse  of  power  P(t).  The  listing  for  the  program  which  has  been  written 
to  numerically  do  this  is  shown  in  Fig.  35a-d.  The  values  of  Eq.  (73)  and  (74)  are 
stored  at  each  mesh  point  (i.  e.  ,  position  z^)  at  a  given  time  t  and  then  updated  to 
time  t  +  dt.  The  beam  tracing  equation  is  then  integrated  forward  starting  from 
z  =  0  into  the  interior  of  the  solid.  Since  we  are  not  able  to  take  advantage  of  the 
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C 

C 


OTrTIME  SA  HP  LE  *NYE  Tr  NO  OF  FRAMES 

™OIm£E?5!"  CALCULATES  the  SELF  FOCUSING  OF  ,  E  ASER  BE«H 

rr::;^01*  •Po-12' 

COHHON/a«NK/  E;?EN“"  ,NETI  ’  "  *-L  .FOOBA  .FILLO. 

szssy&ursiisivzi  •‘5,"ai1  •uR,toi> 

READ  !5  *PARAM  ) 


C 

C 


FOO 

100 


WRITE16.  PARAM) 

namelist/contro/ 

READ  tS.CONTRO? 
WRITE!  6.C0NTR0) 
REA015.100I  G » 
FORHAF  15E10.4) 


OTt NVET.LEGUP  «N  A*  NB  .LETS.  NET.  INTO 


WL.R.AL  . 


AO 


POWER  IN  MEGAWATTS*  COEFFICIENTS  IN  CM/STUT  uuirurTu 

,n,”;/be-;h  »  «•"»  of  foih,^bs„"^’vot;“^^g;: 

EDr£ XP  (-  x ) 

WILL  Yr  (W  L/G.  28  3I8J..2 
plLL  =4  ./  .003»!  l.-EO) 

POO0A=4./.OO3MI  .-tl.-EOJ/XJ 
'PILL0X:4./.(J03*(  1.-E0*!1..X1  )/x 

NE  T I rNET ♦ 1 
GlMr  G 


INCH. RAO 
CM 


OF  CUR 


T=0.0 

IT  IMrN  VET*0T 

WRITEtG.903)  ITIM  .RELAX 
903*r£vM!TrlH1,25x,,TIME  sequences  FOR 

*S* T^m10"  Tl"E  15 

Z(  I  >  rH .FLOAT  ( I -1  j  ♦  0.1E-06 
UlIJrQ.O 


AS  13 


A  !  I )  “AO 
UR  ( I  Jr  0.  o 
5  P(Il:O.Q 
GNEWrO.O 

00  190  LOCKr 1«  NYET 

T=T.OT 

G-GNEW 


nanosecond 


CALL  P  Otf  ER !T  *PAT ) 

GNEWrG IM.PAT 
A(  1 )  rAO 
°(1J=-A0/R 

CAPrUItLY  -8  .*  GNEW 

HEX=  G4.  *E4*  !G  NE  W/ .003  )  *  •? 
MI Xr  35 


/  0.  00  3  *E  2 


00  533  J:l*NETl 
UR!  JjrUR  (J)*EO 
53  3  U!  J) rU !J )*E0 
CALL  YEWIO) 

00  515  Jri.NEl 
IF  !  MIX  .NE.  35  )  GO  TO  832 
IF  f  A  fj  ).  LT.0.0  )  MIX=J 
IF!!P!J)/rj  .g  T.  0.0  j  MIXrj 


PULSES 


Fig.  35a 


Computer  Program 
Equations  for  Beam 


to  Integrate  the  Coupled  Relaxation 
Tracing  with  Finite  Response  Time. 
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33  2  CONTINUE 
ZNr7 I J  3 
UNrUIJ) 

UN1=U(  1 3 
URN- UR  ( J  3 
URNlrURt  J*  1 3 
AN  =  A  (J  I 
PNrPlJ) 

CALL  RAT  El  AN  .PN.  UN  .UNI  »URN  .URN  1.  ZN  3 

A  (  J* 1 3  lA  N 
51  5  PI  J*- 13-  PN 
CALL  Y  EW  1 1  3 

IF  I  HOD  (L  OCK.  LE  TS  3  •  NE  .0  3  GO  TO  190 
WRITE!  6.  2101  PCRIT.GNEW.E2.  E  A.  WL  .R  .AL.  AO  .NET 
210  FORMAT!  35X.  16HCRITICAL  POWER  r.  El  2.  4  •  2HMW /3 7X  .  14HACTUAL  POWER  =• 

2  FJ2.4.2HMW/ 

3  10X  vl  OX  •  SHN2/  N0.t0X.5HN4/N0.5X. 

4  10HWAVELENGTH.  6X» 9HCURVATURE. 

5  9  X  .  6HEX  TEN  T  •  6X.9HRA0IUS  A0.I1X. 

6  4HKESH.  /I  OX  ,6E1  5.  5.  115/  3 
LZ=0 

00  413  MrNA.NB  .LEGUP 
LZ=LZ»1 

413  PO  W !  LZ  3-  GNEW/  3.  14159/AIM)  *»2 

WRITEIG.  200)  V.  IZIN3.NrNA.NB  .LEGU**).  (A  IN  3  •  Nr  NA  *NB  .LEGUP). 

2  IP  IN  3  .  N-NA  .NB  sLESu'M*  I  POW I  N)  .Nr  1  .1 1 ) 

20  0  FORMAT I26X.  ’BEAM  GEOMETRY  WHEN  T=  •  *F 6.  2*  ’NANOSECONDS  * 

J,  /4X8HP0SITI0N.  11 F 1  0.  7/ 6X  .6HR  AO  IU 

3  5.  11  El  0.4/ 7X  .5 HS LOPE  .1  IE  10. 4/4  X.  8HP0YNTING  .1  IE  10  .4  /4X.4HFLUX  3 
LZ -M IX  -6 

00  199  Mrl  .  1 1 
LZ=LZ*  1 

199  POW  I  M)  r  1./3.14159/AILZ)  **2  •  GNEW 

LZ  -  MIX-5 
LZ1=MI  X«-r> 

WRITE!  6.  205)  17  IN  3  •  N=LZ  .LZ1  )•  IA  IN  )•  N=LZ  .LZ1  )• 

2  IP  IN  3 » N-L2  .L  Z1  3.  IP0WIN3.Nrl.il) 

20  5  FORMAT  I26X .  ’GEOMETRY  NEAR  THE  BEAM  WAIST* 

2  /4  X8  HP  OS  IT  ION.  11 T 1  0.  7/  6X  .6HR  AO  IU 

3S.  11  El  0.  4  fit  .5HSL0PE  .1  IE  10.4  /4  X.  8HP0  YNTING  .1  IE  iO  .4  /4X  .  4HFLUX  3 
OPlrPl  MI  «)  -PIMIX-1  3 
ZIPrZI  MI  X)  -P1MIX)  /  OP  1  •  H 

AZIPrA  JMIXNPI  MIX)*!  7IP-ZIMIX)  340«M  /2..IZIP-  Z!MIX))**2/H 
DIGBYr  GNEW/3..14159/AZIP/  AZIP 
WRITE! 6. 899*  Z IP • AZI P. 01 GB Y 

899  FORMAT  I40X.*  WA  1ST  IS  A  T  •  .F 10  .5  .*  CM  •/  l  6X  .’WAIST  RADIUS  IS  ’.F10.5. 

2’CM  .  POWER  FLUX  I  S  •  .E  1 1  .3  .  •  MW/S  0.  CM’///) 

190  CONTINUE 

IN  TO  r  I  NT  0- 1 

IF  I  INT0.E0.1 1  RELAXr5.0 
IF  1 1  NT  0.  GT  .0  3  GO  TO  600 
CALL  EXIT 
END 

Fig.  35b  Computer  Program  to  Integrate  the  Coupled  Relaxation 
Equations  for  Beam  Tracing  with  Finite  Response  Time. 
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PBN-514 


SS!!?™  ??  AT  EC  AN  »PN,  UN  , UNI  »URN  tURN  l.ZN  ) 

COMMON  HEX, CAP  ,H  ,  *L  , *ILL Y 
COMMON /B  AUK/  E2*E<» 

COMMON/BANO/  G»GNEW*X  ,EO  ,NETl  .PTLL.POnp*  DTIIAU 
0  C  U  0  I  —  ||  •  ♦UC>«(E2-E<i*U0))  Mr?.2.,E;*uo»  L 

- 0 (  UO  )*  Y*  UOR 

LETGO:  1 
NOBH 
0(J  AO  rH 

GRAND-  CUN1  -UNI  /OllAO 
GR  IM:(  'JRNi  -URN  )✓ OUAO 
ZEB:ZN 

Ji  5  CONTINUE 
Tl  1 :  H*  PN 

UO-UN,  GR  AND  •  (  ZEB- ZN  ) 

UDR:URN,GRIM*(  ZEB-ZN  ) 

BVrAN 

Cl  INK: POOB  A-GNEV /P Y*  *2 

UO:UO*CLINK 

U0R:U0R*2. -CLINK  /BY--2 
T12=FC  AN  ,U0,  UOR)  -H 
T21  =  H*  (PN-O.  5*  Tl  21 
7EB=ZEB*H/2. 

UO:UN-GRANO  •  (  ZEB-ZN) 

UOR:  URN- GRIMM  ZEB-ZN) 

BY :AN*  0. 5*  Tl  1 

CLINKr  P0  0BA-GNEW/BV--2 

UO:UO- CLINK 

U0R:UDR-2.  -CLINK  /BY--2 
T22:F( BY.UO.UOR) -H 

IFCC  (1  .-T12/T22)  — 2)  .LT.  .005  )  GO  TO  15 
IF((T2  2/PN)**2.LT.0.1E-07)  GO  TO  15 
IFCLETG0.GT.25G.  ANO. NOB.  EO.l  ) 

2 WRITE(6»139)  LETGO  ,  AN 
IF  (NOB  .EO. -1  )  LAMBr-l 
13  3  FORMAT  Cl  X.30H  SLIPPAGE 

2  E15.B)  *  1WLETG0  EOUALS 

IF  (NOB  .EO. -1  )  GO  TO  15 
H:H/2. 

ZEBrZN 

LETGO:  2*  LETGO 
GO  TO  115 
1  5  CONTINUE 
N0B:-1 

T  3 1 :  H-  (P  N»  2.  •  T  22  -T  1 2  ) 

?EB:ZEB*H/2. 

UO:UN»GRANO  •(  ZEB-ZN) 

UDR:UR  N- GRIM-C  ZEB-ZN) 

BY:AN+  2.  -T21  -T  H 
CLINK: POOB  A- GNEV /B  Y- -2 
UO:UO-CL  INK 

U0R:U0R-2. -CLINK  /BY--2 
T32:F( BY tUO.UOR) -H 
ANrAN*  (T  11  ♦«.*  T21*  T31)/6. 

PN:PN-  <  T12,i».  -T22-  T52)/6. 

LETGO:  LETGO- 1 

IF  (LET  GO  .GT.  0)  GO  TO  115 

H:QUAD 

Fig.  35c  Computer  Program  to  Integrate  the  Coupled  Relaxation 
Equations  for  Beam  Tracing  with  Finite  Response  Time. 
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PBN-515 


SUBROUTINE  POWER  IT#  PAT  I 
DIMENSION  TA  Ul  10  I  .  F I  T{  10  I 

DATA  I  TAUI  Nl  #Nrl  ,7 I/O.  .5  .#  7.  #1  0.  #1  5.  #2  0.  #2  5.  /. 

2  I  FIT!  Nl  .N:|  .?  I/O.  *.  7#  1.  #0.8  #0  .5  #0  .1  ./  #M  /7  / 

N-l 

20  N-N^  1 

IFJN.GT.MI  GO  TO  300 
TFITAU  INI. IT  ,T  I  GO  TO  20 

PAT:  FIT!  N-I I  ♦  !  FIT  INI  -FIT  I  N-l  I)  •!  T-TAUI  N-  II  1/  IT  AU  IN  l-TAUIN-  1)  I 

GO  TO  400 
300  «>AT:0.0 
40  0  RETURN 
ENO 


SUBROUTINE  YEWINOKI 

COMMON  /SETL'D/  Al  40  1 1  #U  14  01  I  .AST  14011  #  UR  14  01  1 
COMMON/BAND/  G.GNEW.X  #E0  .NET  1  .  PILL  .POOBA  .P  ILLOW 

IFlNOK.EO.il  GO  TO  99 
00  20  J=1.NET1 

UR  1  J  ):  UR  I J  I  ♦  2  .  *P  ILLOW*  G  /A  IJ  1*  *4 

20  U<  Jl  :U  IJ  1»PILL0V*G/  AIJI**2 
RETURN 

99  00  35  J=1.NET1 

URIJ1:URIJI»2.  *POOBA  *  GNEW /A  IJ  !•  *4 
35  U I  J I  :U  I J I  ♦  POOBA  •  GNEW/AIJ}**2 

RETURN 
ENO 


Fig.  35d  Computer  Program  to  Integrate  the  Coupled 

Relaxation  Equations  for  Beam  Tracing  with 
Finite  Response  Time. 
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fictitious  potential  approach  in  this  program,  some  experimentation  is  usually 
necessary  in  choosing  the  mesh  size  in  the  z-direction  small  enough  (the  mesh 

size  is  H  =  AL/NET  where  AL  is  the  slab  thickness)  in  order  to  prevent 
instabilities  of  the  "giant  planet"  type. 


Some  results  are  shown  in  Sec.  IV -G  comparing  a  solid  with  a  response 
time  of  5  nanoseconds  and  one  having  an  instantaneous  response  time. 


G.  Numerical  Examples 


The  program  was  carried  out  for 
parameters 


some  examples  using  the  following 


ho  ■  1-0 

h2  =  2. 0  x  10  cm2/statvolt 

h4  =  3.2  X10  24  cm4/  statvolt2 

where 

h  =  h0  +  h2u  -  h4u2 

Figure  36a  shows  the  pulse  shape,  which  approximates  an  experimental 
pulse  while  Fig.  36b  shows  the  focal  position  as  a  function  of  time  for  an  entrance 
face  beam  radius  of  0.184  mm  and  a  beam  convergence  angle  of  73.  6  milliradians 
which  would  imply  a  geometrical  optics  focus  at  0.25cm.  This  corresponds  to  an 
experiments1  geometry.  The  relaxation  time  tq  was  taken  to  be  a  very  small 
GO  seconds)  compared  to  the  pulse  time. 

In  Fig.  37  we  show  the  computed  focal  position  for  tq  =  5  nanoseconds 
for  the  same  pulse  shape  and  physical  parameters  used  in  Fig.  36.  In  this 
example  the  entrance  face  beam  radius  was  taken  as  0.05mm  and  the  convergence 
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Position  of  focal  point  -  cm  Power  (arbitrary  scale) 


PBN-503 


Fig.  36  Calculated  Distance  of  the  Focal  Point  From  Entrance  Face 

as  a  Function  of  Time  for  an  Initial  Beam  Radius  of  0. 184  mm 
and  a  Convergence  Angle  of  73.6  Milliradians. 
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Position  of  foccl  point  (cm)  Power  (arbitrary  scale) 


PBN-504 


Fig.  37  Calculated  Distance  of  the  Focal  Point  from  Entrance  Face  as 
a  Function  of  Time  for  an  Initial  Beam  Radius  of  50  urn  and  a 
Convergence  Angle  of  20  Milliradians. 
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angle  as  20  milliradians.  The  focal  length  moves  in  toward  the  entrance  surface 
as  the  power  increases  but  the  shortest  focal  length  is  attained  considerably  after 
the  power  reaches  its  maximum  in  contrast  to  the  result  for  an  instantaneous 
response. 


The  beam  radius  as  a  function  of  position  is  plotted  in  Fig.  38  for  the 
case  of  instantaneous  relaxation  at  a  time  4  nanoseconds  after  the  start  of  the 
pulse  when  the  instantaneous  toSal  power  is  1.18MW.  Since  P  is  greater  than 
Pcrit  (where  pcrit  is  deftned  as  the  power  at  which  the  quadratic  term  in  E2  is 
balanced  by  the  diffraction  term),  the  beam  is  initially  self-foe  ’sed,  but  when 
beam  radius  is  less  than  3  microns  the  r?4  term  finally  dominates  and  causes 
the  beam  to  diverge  once  more.  Trapping  is  much  less  effective  for  the  form 
given  by  Eq.  (28)  than  it  is  for  the  more  realistic  index  of  refraction  given 
by  Eq.  (29),  which  gives  the  correct  saturation  limit. 

Another  situation  of  experimental  interest  was  a  collimated  beam 
which  had  an  entrance  beam  radius  of  313pm.  With  a  power  of  1. 97  MW,  which 
was  nearly  four  times  the  critical  power  of  0.  53  MW,  the  convergence  was 
extremely  slow  and  the  beam  had  converged  only  to  300pm  after  10  cm  travel, 
a  beam  length  large  compared  with  the  actual  sample  size.  This  result  was 
in  agreement  with  the  experimental  results,  in  that  no  indication  of  self-focusing 
could  be  found  with  this  beam  geometry. 


108 


U8  mW  Beam 


(^)  srijpo-l  UJD8g 


109 


V. 


CONCLUSIONS 


The  objective  of  this  program  is  to  determine  how  and  why 
permanent  damage  is  produced  and  to  find  means  vo  minimize  or  avoid 
its  deleterious  effects  on  laser  systems.  The  work  to  date  can  be 
grouped  into  three  general  areas:  1,  studies  of  the  residual  damage 
after  irradiation  using  optical  and  electron  microscopy,  2,  studies 
of  the  dynamics  of  the  damage  process  using  high  speed  streak 
photography  and  3,  determination  of  the  probability  that  a  particular 
material.  Th€'  following  list  summarizes  the  major  results  reached 
to  date: 


1.  Residual  damage  was  found  only  after  the  laser  irradiation 
caused  a  visible  spark  when  passing  through  the  sample. 

2.  The  importance  of  uniform  ("hot  spot"  free)  illumination  in 
increasing  the  ability  of  a  material  to  withstand  laser  irradiation  was 
demonstrated  by  comparing  TEMqq  and  multimode  damage. 

3.  Contrary  to  previously  held  opinions,  when  LiNbOg  and  KDP 
are  used  as  Pockel' s  cell  or  frequency  doubler  crystals,  internal  damage 
is  produced  at  the  same  or  lower  levels  of  irradiation  than  those  required 
to  produce  surface  damage. 

4.  The  laser  beam' s  focusing  conditions  determine  the  form  of 
the  residual  damage;  that  is,  whether  it  is  internal,  entrance  or  exit 
face  damage. 

5.  Internal  filamentary  damage  in  LiNbOg  was  found  to  be  composed 
of  a  series  of  very  fine,  nearly  planar  cracks  which  intersect  to  form  one 
or  more  long  lines  having  diameter  =  0.  4|im.  When  illuminated  with 
visible  light,  diffraction  effects  make  these  lines  appear  blurred  and  give 
the  impression  of  a  "  tube"  of  damaged  material. 
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6.  Expei  iments  designed  to  study  the  dynamics  of  damage 
formation  by  photographing  the  laser  induced  breakdown  (the  spark) 
with  an  image  converter  streak  camera  were  initiated.  The  results 
of  these  measurements  showed  that  when  light  is  focused  inside  the 
material  the  first  damage  to  occur  is  near  the  focal  point  and  additional 
damage  occurs  upstream  at  later  times.  All  laser-induced  internal 
breakdowns  were  initiated  during  intervals  of  time  when  the  level  of 
irradiation  increased  or  was  maximum,  and  never  when  the  intensity 
decreased. 

7.  A  model  for  the  results  in  6  was  developed  in  which  self 
focusing  was  identified  as  likely  to  be  responsible  for  the  upstream 
movement  of  the  sparks  and  a  fast  response  self-focusing  process  was 
shown  to  be  able  to  explain  the  fact  that  internal  breakdowns  were  only 
produced  when  the  pulse  intensity  increased.  The  importance  of  beam 
geometry  and  non-self-focusing  configurations  were  investigated. 

8.  There  is,  in  general,  no  sharp  threshold  for  laser  induced 
damage.  There  is,  however,  some  probability  that  a  particular  pulse 
will  damage  a  particular  material.  At  very  high  levels  of  irradiation 
this  probability  is  unity,  at  lower  levels  it  is  less  than  one  and  at  very 
low  levels  the  probability  for  damage  may  be  vanishingly  small.  This 
was  studied  with  both  Nd:YAG  and  ruby  lasers. 

9.  Measurements  of  the  dependence  of  damage  probability  on  the 
laser  beam' s  power  density  for  ten  different  materials  have  led  to  a 
model  xor  the  damage  mechanism  based  on  avalanche  breakdown.  Good 
qualitative  agreement  has  been  obtained. 

10.  Using  the  streak  camera,  preliminary  measurements  of  the 
distribution  of  breakdown  starting  times  as  a  function  of  damage  probability 
have  been  obtained.  These  data  show  that  the  most  likely  time  for  a  break¬ 
down  to  occur  is  before  the  peak  of  the  laser  pulse.  This  data  may  be 
explained  using  the  probabilistic  interpretation  of  laser  induced  damage. 
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APPENDIX  A 


THERMAL  SELF-FOCUSING  IN  LiNbO^* 

The  self-focusing  mechanisms  for  solids  most  often  mentioned  are 
2 

1)  electro -striction  and  2)  the  hermal-heating-induced  index-of-refraction 
4 

changes.  In  the  former,  the  index  of  refraction  in  the  irradiated  region 
increases  due  to  the  electrostrictively  induced  stress  and  the  beam  is 
focused.  In  the  latter  if  the  index  of  refraction  increases  with  temperature, 
and  if  there  is  some  absorption  process  whereby  the  medium' s  temperature 
can  be  increased,  then  the  beam  can  be  focused. 

v. 

To  see  if  self-focusing  due  to  thermal  effects  can  be  significant  in 
LiNbOg,  the  absorption  in  this  material  at  the  laser  wavelength  had  to  be 
measured.  If  the  absorption  per  centimeter  in  the  material  is  |3 ,  the 
volume  V  will  be  heated  by  AT  =  &£/VC  by  the  passage  of  a  pulse  of 
light  of  energy  £.  C  is  the  specific  heat  in  joul  c  per  cm3  per  °K  and 
the  pulse  duration  is  assumed  so  short  that  thermal  conductivity  can  be 
neglected. 

In  LiNbOg  the  focal  volume  of  a  typical  lens  used  in  this  work 
will  be  heated  one  centigrade  degree  by  a  1  mJ  pulse  if 

/3  =  0.001  cm  *  . 

Since  so  small  an  absorption  cannot  be  measured  using  conventional 
spectrophotometry  (i.e. ,  in  a  device  such  as  a  Cary  14  Dual  Beam 
Spectrophotometer)  the  following  very  sensitive  procedure  was  developed. 

If  heating  due  to  absorption  at  1.06  fim  occurred,  than  a  LiNbOg  crystal 
which  was  used  as  an  "a-axis"  frequency-doubler  would  require 
continuous  temperature  adjustment.  If  such  a  crystal  were  fixed  at  a 


The  measurements  described  in  this  section  v/ere  performed  as  part 
of  work  for  the  Electronics  Research  Center  under  NASA  Contract 
NAS  12-21 55. 
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temperature  below  the  index -matching  temperature  T  and  then  exposed 
to  light  at  1.06  pm,  the  heating  due  to  absorption  could  be  detected  by 
monitoring  the  increase  in  harmonic  output  as  the  temperature  of  the 
irradiated  volume  rose  to  Tm.  Since  the  initial  temperature  and  the 
index-matching  temperature  are  known,  the  absorption  can  be  estimated 
by  measuring  the  incident  laser  power  and  the  time  required  to  reach  T  . 
This  type  of  experiment  was  performed  using  "hot"  LiNbOg  and  the 
resulting  data,  the  expected  "sin  Kt/Kt",  function  is  shown  in  Fig.  A-l. 

A  10-pps  pulsed  Nd:YAG  laser  was  used.  The  average  power  was  1W 
at  1.06  pm,  the  maximum  harmonic  power  was  0.005  W  and  the  irradiated 
volume  was  a  cylinder  1.5  cm  long  and  ~  0. 5  cm  in  diameter.  When  the 
average  laser  power  was  ~  0, 1  W  and  the  harmonic  power  was  ~  0.02  W, 
no  appreciable  heating  was  observed.  Thus  the  heating  is  due  to  1.  06  pm 
absorption,  and  so  by  neglecting  the  thermal  conductivity  of  the  material 
t  :  smallest  absorption  which  could  give  the  requisite  temperature  increase 
is  found  to  be  |3  *  0.002  cm"1. 

The  possibility  that  thermal  effects  can  lead  to  self-focusing  in 
LiNbOg  crystals  therefore  cannot  be  ruled  out.  In  fact,  since  most 
application  involve  laser  pulses  o  ‘  many  mili joules,  this  possibility 
must  be  included  in  any  analysis  o  ?  damage  in  LiNbOg.  Similar 
consideration  must  be  given  other  materials  of  interest. 

It  is  appropriate  to  ask  whether  the  absorption  described  above 
should  account  for  the  melting  observed  in  Figs.  12-14.  The  melting 
point  of  LiNbOg  is  1260*C  and  the  absorption  coefficient  at  1.06  pm  is 
“  0.002  cm  .  Thus  the  entire  pulse  energy  would  have  to  be  absorbed 
in  the  0.006  cm  diameter  by  ~  1  pm  deep  crater  in  Fig.  12,  to  achieve 
a  temperature  increase  to  the  melting  point  of  LiNbOg.  However,  since 
the  absorption  coefficient  is  so  small,  it  is  extremely  unlikely  that  the 
pulse  can  be  completely  absorbed  in  the  first  1  pm  of  material. 
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